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The formulation and resolution of integrable lattice statistical models in a 
quantum group covariant way is the subject of this review. The Bethe Ansatz 
turns to be remarkably useful to implement quantum group symmetries and to 
provide quantum group representations even when q is a root of unity. We start by 
solving the six-vertex model with fixed boundary conditions (FBC) that guarantee 
exact SU(2) q invariance on the lattice. The algebra of the Yang-Baxter (YB) and 
SU(2) q generators turns to close and the transfer matrix is SU(2) q invariant for 
FBC. In addition, the infinite spectral parameter limit of the YB generators yields 
cleanly the SU(2) q generators. The Bethe Ansatz states constructed for FBC arc 
shown to be highest weights of SU(2) q . The light-cone evolution operator for 
FBC is introduced and shown to follow from the row-to-row FBC transfer matrix 
with alternating inhomogeneities. This operator is shown to describe the SOS 
model after an appropiate gauge choice. Using this FBC light-cone approach, the 
scaling limit of both six-vertex and SOS models easily follows. Finally, the higher 
level Bethe Ansatz equations (describing the physical excitations) are explicitly 
derived for FBC. We then solve the RSOS(p) models on the light-cone lattice with 
fixed boundary conditions by disentangling the type II representations of SU(2) q , 
at q = e l7V / p , from the full SOS spectrum obtained through Algebraic Bethe Ansatz. 
The rule which realizes the quantum group reduction to the RSOS states is that 
there must not be singular roots in the solutions of the Bethe Ansatz equations 
describing the states with quantum spin J < (p — l)/2. By studying how this 
rule is active on the particle states, we are able to give a microscopic derivation 
of the lattice S— matrix of the massive kinks. The correspondence between the 
light-cone Six-Vertex model and the Sine-Gordon field theory implies that the 
continuum limit of the RSOS(p+ 1) model is to be identified with the p— restricted 
Sine-Gordon field theory. 
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1. Introduction 



As is by now well known integrability is a consequence of the Yang-Baxter equa- 
tion (YBE) in two-dimensional lattice models and two-dimensional quantum field 
theory (QFT), (for recent reviews see for example ref. [1]) . More precisely, a statisti- 
cal model is integrable when the local weights are solutions of the YBE. Analogously, 
for two dimensional integrable QFT the two-body S-matrix fulfils the YBE. 

Quantum groups are closely related to Yang-Baxter algebras [1,2]. However, quan- 
tum group invariance holds for an integrable lattice model only for specific choices of 
the boundary conditions. As we showed in refs.[3-4] (see also refs. [5,6,7]), choosing 
fixed boundary conditions (FBC), the transfer matrix commutes with the quantum 
group generators. 

The purpose of this paper is to review the work in collaboration with Claudio 
Destri in refs. [3-4] on integrable lattice models and their scaling limit using a fully 
quantum group covariant Bethe Ansatz (BA) framework. Let us recall that the YB 
algebra of monodromy operators acting on the space of physical states, is the main 
tool to construct the transfer matrix eigenvectors by the (algebraic) Bethe Ansatz. To 
do that we choose FBC (Dirichlet type) and use the Sklyanin-Cherednik[8,9] construc- 
tion of the Yang-Baxter algebra. In this framework, besides the R-matrix defining 
the local statistical weights, there are two matrices K ± (9) that define the boundary 
conditions. K ± {9) must fulfil a set of equations [eqs.(2.12) and (2.15)] in order to 
respect integrability. (FBC is one special case out of a continuous family of boundary 
conditions compatible with integrability). 

The appropiate monodromy operators U a b(\,£j) for a N-sites line take now the 
form depicted in fig.l where arbitrary inhomogeneities uj{ , (1 < i < N), are allowed 
at each site. We compute the large 9 limit of the monodromy operators U a i,(\, Co) 
and find that they are just the quantum group generators. We do that explicitly for 
the six- vertex model where the quantum group is SU(2) q . We find in this way an 
explicit representation of the SU (2) q generators acting on the space of states. 

In addition, the 9 — > oo limit of the Yang-Baxter algebra for the U a b(\,uj) shows 
that the transfer matrix t(X,oj) commutes with the SU(2) q generators and that the 
algebra of the U a }>(\, uj) (1 < a,b < 2) with the quantum group generators ( J± and 
q Jz ) closes [see eqs.(2.33)-(2.37)]. Moreover, t(X,uj) in the 9 ^ oo limit yields the 
q-Casimir operator C q through 

t(oo,Cu) = q + q~ 1 + (q-q- 1 ) 2 C q (1.1) 

It must be recalled that for periodic boundary conditions (PBC) the 9 — > oo limit 
of the Yang-Baxter algebra also gives the quantum group generators but that the 
algebra with the PBC monodromy operators does not close [1] . 
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Then, we investigate the BA construction in this quantum group covariant frame- 
work. It is natural to define for FBC a creation operator of pseudoparticles B(9) 
(proportional to U\2{\w) ) which is odd in 9 [eq.(3.3)]. The exact eigenvectors of 
the transfer matrix t(X, cu) are then given by 

= B(v 1 )B(v 2 )....B(v r )n (1.2) 

where v 1, V2, ■ ■ ■ , v r , fulfil the BA equations (3.4) (BAE) and Q is the ferromagnetic 
ground state (3.2). We show that only BAE roots with strictly positive real part 
must be consider (Revj > 0). In particular, vj = and purely imaginary roots 
must be discarded. We map the BAE (3.4) for FBC onto BAE for PBC in 2N sites 
[see eq.(5.2)]. We find the usual PBC BAE plus an extra source and two important 
constraints : 

a) the total number of roots is even and they are symmetrically distributed with 
respect to the origin, 

b) a root at the origin as well as purely imaginary roots are excluded. 

Therefore, the antiferroelectric ground state (and the excitations on the top of 
it) contain always a hole at the origin. This hole combined with the extra source 
accounts for the surface energy (see for example ref.[12]). 

Starting from the general BA state (1.2) we prove that they are highest weights 
for SU{2) q . That is, 

J + ^(v) = (1.3) 

provided the BAE (3.4) hold for i>i,i>2, . . . ,v r . Therefore, the eigenvalues of t(6,u) 
are degenerate with respect to the quantum group and the eigenvectors : 

J^(v), (J-) 2 ^(v),..., (J-) 2J ^(v), (1.4) 

are linearly independent from ty(v) . 

SOS and vertex models are related by the vertex-face correspondence. (The 
degrees of freedom lie on faces for SOS models whereas they lie on links for the 
vertex model). The correspondence between them amounts to an application of the 
g-analog of the Wigner-Eckart theorem. As we explain in sec. 4, the SOS space of 
states is identical to the set of maximal weight six vertex states in a quantum group 
invariant framework like ours. It then follows from eq.(1.3) that the Solid-On-Solid 
(SOS) states are just the f.b.c. BA states given by eq.(1.2). The six vertex Hilbert 
space includes the whole SU(2) q multiplets and follow then by repeatedly applying 
the lowering operator J_ to the highest weight BA states (1.2) as shown in eq.(1.4). 
In other words, we have found the BA solution of the SOS model since we derived 
the transfer matrix eigenstates. In ref.[10] an alternative but equivalent solution of 
the SOS model is derived using a BA in face language. In addition PSOS (periodic 
SOS models where the face states I and I + p are identified) are solved in ref.[10]. 
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The light-cone approach is a direct way to give a field theoretical interpretation to 
a lattice model and furthermore obtain its scaling limit as a massive QFT. The light- 
cone approach with periodic boundary conditions has been investigated in refs. [19-20] 
. Here we consider this approach in the case of fixed boundary conditions leading to 
a quantum group invariant framework. 

In this context we show in sec. 4 that the diagonal-to-diagonal transfer matrix 
U (©) [fig. 6] can be obtained from the row-to-row transfer matrix t(X, Co) choosing the 
inhomogeneities appropiately. This is analogous to the relation found in ref. [19,20] 
for PBC. From the transfer matrix U (©) we define the lattice hamiltonian through 

if = - log 17(6) (1.5) 

CL 

where a is the lattice spacing. In the a — > 0, — > oo limit this operator defines the 
continuous QFT hamiltonian provided the renormalized mass scale [see eq.(4.17)] is 
kept fixed. The evolution operator U (©) can be considered both in the vertex or in 
the face language. In face language, it has a simple expression provided we make an 
appropiate gauge transformation. That is, if one transforms each local i?-matrix to a 
matrix R such that the SU(2) q symmetry holds locally. In this way, we show that the 
associated light-cone evolution operator U(Q) just describes the SOS ABF model[13] 
[eq.(4.27)]. 

As it is known, the physical states above the AF vacuum are described by the 
higher level BAE[14]. We obtain the higher level BAE for FBC in sec. 5 [eq.(5.14)]. 
They are the starting point to study the excitations in SOS and RSOS models . To 
conclude we show that the BAE for FBC admit solutions at infinity only when 7/V 
is a rational number. This is precisely the case when RSOS models can be defined 
and when the representations of SU(2) q algebras cease to be isomorphic to usual 
SU(2). Recall that when 7/V is rational, type I representations are reducible but 
indecomposable, type II are irreducible as in SU(2) (see for example ref. [5,6]). 

It is known that the Six- Vertex (6V) model, in the so-called light-cone formula- 
tion and with periodic boundary conditions (p.b.c), yields the Sine-Gordon massive 
field theory in an appropriate scaling limit [19]. Hence the light-cone 6V model can 
be regarded as an exactly integrable lattice (minkowskian) regularization of the SG 
model. 

Recently, the hidden invariance of the SG model under the quantum group SU(2) q 
was exhibited [21]. Our quantum group invariant light-cone formulation, provides a 
lattice formulation where such hidden invariance appears starting from first principles. 
That is, not on a bootstrap framework but deriving the field theory as a rigorous 
scaling limit of the six-vertex model. 

We present in section 6 and 7 the derivation of the factorized S— matrices on 
the lattice, i.e. still in the presence of the UV cutoff. This derivation is based on 
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the "renormalization" of the BA Equations, which consists in removing the infinitely 
many roots describing the ground state. What is left is once again a f.b.c. BA 
structure involving the lattice rapidities of the physical excitations (the particles of 
the model) and the roots of the higher-level BAE obtained in paper I. The explicit 
form of two-body 5 1 — matrix for the 6V model and the SOS model can be extracted 
in a precise way (eft. eq. (6.11)) from this higher-level BA structure. In the massive 
scaling limit these lattice scattering amplitudes become the relativistic S— matrices 
of the SG model (or Massive Thirring model) and of the continuum SOS model. 
Let us remark that the SOS S— matrix, although closely related to the 6V and SG 
S— matrices from the analytical point of view, is conceptually different. It describes 
the scattering of kinks interpolating between renormalized local vacua labelled by 
integers. This kink S— matrix is most conveniently expressed in the so-called face 
language (see eq.(7.3)). 

In section 8 we investigate the f.b.c. BAE (eq.(2.4)) when the quantum group 
deformation parameter q is a root of unity, say q p = ±1, with p some integer larger 
than 2 (the case p = 2 being trivial). In this case it is known that RSOS(p) model 
can be introduced by restricting to the finite set (1, 2, ... ,p — 1) the allowed values 
of the local height variables of the SOS model [13]. At the same time when q p = ±1, 
the representations of SU(2) q algebras cease to be isomorphic to usual SU(2). Recall 
that when •y/n is rational, type I representations are reducible but indecomposable, 
type II are irreducible as in SU(2) (see for example ref.[5,6]). The restriction leading 
to the RSOS model from the SOS model is equivalent to the projection of the full 
SOS Hilbert space (which is formed by the highest weight states of SU(2) q ) to the 
subspace spanned by the type II representations [5,6]. That is, those representations 
which remain irreducible when q becomes a root of unity. Our results on this matters, 
in the BA context, can be summarized as follows: 

a) Only when q is a root of unity, the f.b.c. BAE admit singular roots (that is 
vanishing z— roots or diverging t>— roots, in the notations of sec.l). 

b) When q tend to a root of unity, say q p = ±1, the BAE solutions can be divided 
into regular and singular solutions, having, respectively, no singular roots or 
some singular roots. Regular solutions correspond to irreducible type II repre- 
sentations. Singular solutions with r singular roots correspond to the reducible 
and generally indecomposable type I representations obtained by mixing two 
standard SU(2) q irreps of spin J and J + r (we recall that J = N/2 — M, 
where N is the spatial size of the lattice and M the number of BA roots). Then 
necessarily r < p. 

c) The r singular roots z\, Z2, ■ ■ ■ , z r vanish with fixed ratios 

z - = J~ 1 Zl , uj = e 2ni / r , \<j<r (1.6) 

In terms of the more traditional hyperbolic parametrization, with v— roots re- 
lated by Vj = — ilogZj to the z— roots, the singular roots form an asymptotic 
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string-like configuration. They have a common diverging real part and are 
separated by ir/r in the imaginary direction. 

So we see that the RSOS eigenstates are easily singled out from the full set of 
BA eigenstates of the 6V or SOS transfer matrix. One must retain all and only 
those BAE solutions with M > (N — p + l)/2 (which correspond to states with 
J < (p — l)/2) and with all M z— roots different from zero. This provides 
therefore an exact, explicit and quite simple solution for the RSOS model on the 
lattice (with suitable boundary conditions, as we shall later see). In particular the 
ground state of the 6V, SOS and RSOS models is the same f.b.c. BA state (in the 
thermodynamic limit N — > oo at fixed lattice spacing). It is the unique SU(2) q singlet 
with all real positive roots and no holes. The local height configuration which, loosely 
speaking, dominates this ground state can be depicted as a sequence of bare kinks 
jumping back and forth between neighboring bare vacua (see fig. 10). In the massive 
scaling limit proper of the light-cone approach [19], this ground state becomes the 
physical vacuum of the SG model as well as of and all the Restricted SG field theories. 

In the discussion ending sec. 8, we argue that the kink S— matrix for the excita- 
tions of the RSOS models follows indeed by restriction from that of the SOS model. 
In the scaling limit it is to be identified with the relativistic S— matrix of the Re- 
stricted SG models. All these field-theoretical S— matrices are naturally related to 
the Boltzmann weights of the respective lattice models. This is because the higher- 
level BAE are identical in form to the "bare" BAE, apart from the renormalization 
of the anisotropy parameter 7 (related to q by q = e 17 ) 

7-7= — (1-7) 

7T — 7 

and the replacement of the rapidity cutoff ±0 with the suitably scaled rapidities 9j 
of the physical excitations 

±©^^L_ (1.8) 

7T — 7 

In particular, since 7 = n/p for the RSOS(p) model, eq.(1.7) yields p — > p — 1. This 
shows that the renormalized local vacua £ run from 1 to p — 2 when the bare local 
heights £ n run from 1 to p — 1. The higher-level BA structure of the light-cone 6V 
model (or lattice regularized SG model) thus provides a microscopic derivation of the 
bootstrap construction of ref. [21] and explains why the S— matrix of the p— restricted 
SG field theory has the same functional form of the Boltzmann weights of the lattice 
RSOS(p+ 1) model. Moreover, the well-known correspondence between the critical 
RSOS(p) models and the Minimal CFT Models M p imply the natural identification of 
the massive p— restricted SG model with a completely massive relevant perturbation 
of M p . This is generally recognized as the perturbation induced by the primary 
operator 0^3 with negative coupling. 
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Two detailed examples of the BA realization of the quantun group reduction to 
the RSOS models are presented in section 9. We considered the simplest cases p = 3 
and p — 4. The RSOS (3) model is a trivial statistical system with only one state, since 
all the local heigths are fixed once we choose, for example, the boundary condition 
£q = 1. In our quantum group covariant f.b.c. construction this corresponds to the 
existence of one and only one type II state when 7 = n/3. We then obtain the 
following purely mathematical result: for any N > 2 and real w = exp(— 20) the set 
ofBAE 

ZjW - e™/ 3 Z j - we ni / 3 \ N _ [ S ] Zj - z k e 2ni / 3 Zj z k - e 2ni / 3 < . < 

Zjwe™/ 3 - 1 zje™/ 3 -w j ' zje 2 ™/ 3 - z k Zj z k e 2ni / 3 - 1 ' ~ 3 ~ 

(1-9) 

admit one and only one solution with non-zero roots within the unit disk \z\ < 1. In 
addition, these roots are all real and positive. 

The RSOS (4) model can be exactly mapped into an anisotropic Ising model as 
showed in eqs. (9.1-5). In our case the horizontal and vertical Ising couplings turn 
out to be O— dependent complex numbers. For even N, the Ising spins are fixed on 
both space boundaries. For odd iV the spins are fixed on the left and free to vary 
on the right. We analyzed the BAE for the RSOS(4) model in some detail. In the 
thermodynamic limit the ground state, as already stated, is common to the 6V and 
SOS moodels. The elementary excitations correspond to the presence of holes in 
the sea of real roots charactering the ground state. Each holes describes a physical 
particle or kink and may be accompanied by complex roots. In sec. 9 we argue that 
in the RSOS (4) case a state with v holes necessarily contains \y/2\ two-strings those 
position is entirely fixed by the holes. Notice that here the number v of holes can be 
odd even for even N. This is not the case for the 6V or SOS models, where v is always 
even for N even. What happens is that when 7 — > (7r/4) _ the largest real v— root 
diverges in the J = 1 states of the 6V and SOS models. Therefore, these states get 
mixed with J = 2 states into type I representations and do not belong to the RSOS (4) 
Hilbert space. The RSOS (4) states with J = 1 are obtained by choosing the largest 
quantum integer In/2-1 = (see the Appendix). There is no root associated to 
N/2 + 1. It follows that these states, from the 6V and SOS viewpoint, have a cutoff 
dependent term in the energy equal to ir/a, where a is the lattice spacing (loosely 
speaking, one could say that "there is a hole at infinity"). They are removed from the 
physical SG spectrum in the continuum limit. We are thus led to propose as RSOS(4) 
hamiltonian, for even N 

HrsoS(4) = Hsosil = tt/4) - a^nJ (1.10) 

where H sos is given by eqs. (4.14), (4.10), (4.21-27) and J = or 1. In this 
way, the particle content of the light-cone RSOS (4) and the corresponding S'-matrix 
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turn out to coincide with the results of the bootstrap-like approach of ref. [21]. 
Eq. (1.10) defines, in the scaling limit, the hamiltonian of the (p = 3)— restricted SG 
model. Notice, in this respect, that the higher-level BAE (5.17) and (6.10) completely 
determine the physical states in terms of renormalized parameters. 

To summarize, the picture we get from the BA solution of the f.b.c. 6V, SOS and 
RSOS(p) lattice models is as follows. Performing the scaling limit whithin the light- 
cone approach, these lattice models yields respectively: the SG model (or Massive 
Thirring Model), a truncated SG and the (p — l)restricted SG models. For the SG 
model we have essentially nothing to add to the existing literature, apart from the 
explicit unveiling, at the regularized microscopic level, of its SU (2) q invariance and for 
a better derivation of the S— matrix. The truncated SG follows by keeping only the 
highest weight states with respect to SU (2) q , that is the kernel of the raising operator 
J + . Finally we showed that the RSOS(p) lattice models with trigonometric weights 
yield in the scaling limit proper of the light-cone approach the (p — l)restricted SG 
field theories formulated at the bootstrap level in ref. [21]. 

2. Boundary conditions in lattice integrable models 

Let us consider an integrable vertex model with R-matrix R^(9) (see fig.l). Each 
element R^(9) defines the statistical weight of this configuration. We assume R{9) 
to fulfil the Yang-Baxter equations. The indices a, b, c, d, are assumed to run from 1 
to q with q > 2. 

[1 <g> R{9 - 6')] \R{6) 1] [1 <g> R{0')] = [R(9') <g> l] [1 <g> R(6)] [R(6 - 9') <g> l] 

(2.1) 

We shall assume T and P invariance for R{9) 

RfM = Ratio) = R b d a M (2.2) 

In addition, we assume R(9) to be regular, that is 

R(0) = cl or R?M = c6 a c 5 b d (2.3) 

where c is a numerical constant. Eqs. (2.1) and (2.3) imply the unitarity-related 
equation 

R{9) R{-9) = p{9) 1 (2.4) 

where p{9) is an even c-number function. Furthermore, we assume crossed-unitarity. 
That is, 

R{9) R(-9 - 2rj) = p{9) 1 (2.5) 
where 77 is a constant and R^{9) = Rl^{9) . 
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Let us consider now a NxN' square lattice with periodic boundary conditions. 
Then the row-to-row transfer matrix is given by 



T{e,u) = Y J Taa{e,u) (2.6) 

a 

where the operators T ab (9, Co) defined by (see fig. 2) 

T ab (9,Co) = W(0-^i)®*o2ai(0-^)® ® W-i(^-^jv) (2.7) 

ai,...,ojv-i 

act on the vertical space V = ®i<i<jv ^ > ^ = ^ > an< ^ the simplest choice for 

the local vertices is [t a b(9)] c d = ^co(^)- m ec L-(2.7) coi,t02, .....,un stand for arbitrary 
inhomogeneity parameters. The T ab (9,Co) fulfil the YB algebra 

i?(A - ji) [T(A, Co) <g> T(n, Co)] = [T(n, Co) ® T(A, a;)] R(X - //) (2.8) 

as follows from eq.(2.1). Thus, the t(9,Co) are a commuting family 

[r(9,Co),r(9 , ,Co)} =0 (2.9) 

Let us now consider the generalization to other boundary conditions compatible with 
integr ability [8]. Define (see fig. 3) 

U ab (9,Co) = Y,T ad {0,Co)K- c {e)T-\-e,Co) (2.10) 

cd 

where T~ 1 (—9, Co) is the inverse in both the horizontal and vertical spaces : 

J2T ab (9,Co)T bc 1 (-9,Co)= 5 a c l (2.11) 
b 

and I is the identity on the vertical space V . Summation over indices of the vertical 
space V are omitted both in eqs.(2.10)and (2.11)(cft. fig. 3). K~{9) in eq.(2.10) is a 
qxq matrix solely acting on the horizontal space. It must fulfil [8] 

R(\-ri[K-(\)®l]R(\ + ri [K-(jm)®1] = 

\K~{h) <S> l] R(X + n) \K~{\) <g) l] R{\ — n) 

(Notice, that our R-matrix differs from ref. [8] in a permutation matrix R — > PR, 
P£(6) = 5«5 b c ). 
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As is well known, R^(6 — 9') has the interpretation of scattering amplitude for a 
two-body collision where a(d) and 6(c) label the initial (final) states of two particles 
with rapidities 9 and 9' respectively. In this S-matrix context, K~ h {9) is the scattering 
amplitude for one particle with a rigid wall on the left, a and b labelling the initial 
and final states and 9 being its final rapidity (see fig.4). 

Thanks to eqs.(2.8) and (2.12), U(9,Co) fulfils the Yang-Baxter algebra 



R(X - (j,) [U(X, Co) <g> 1]R(\ + fx) [U(fx, Co) <g> J] = 

[Ufa, Co) <g> 1] R(\ + fi) [U(X, Co) <g> 1] R(\ - fi) 



(2.13) 



The transfer matrix is given now by 



t(\,Co) = J2K b ^ + v)Uab(\,Co) 

ab 



(2.14) 



Here K + (X) describes the scattering with a rigid wall on the right (eft. fig. 5). It is a 
solution of the equation : 



R(X K+(X)]R{X + K+{li)} = 

[1 ® K+{fi)] R{X + /x) [1 <g) K+{X)] R{X - n) 

It follows from eqs.(2.13)-(2.15) that t(X,Co) is a commuting family 

[t(6,u>), t(9',Co)] = 



(2.15) 



(2.16) 



As we can see from eqs. (2.10) and (2.14) the boundary conditions associated to 
t(9,Co) follow from the form of -^ + (A) and K~(X) sitting on the right and left borders, 
respectively. 

From now on we shall consider the case of the six-vertex model, where q = 2 and 
the R-matrix reads, in terms of a generic spectral parameter 9, 



R{9) 



/I 








o\ 





c 


b 








b 


c 





u 








1/ 



b = b(9, 1 ) = - 
c = c(0, 7 ) 



sinh^ 



sinh(«7 — 9) 
sinh i'-f 



(2.17) 



sinh(«7 — 9) 



where the anisotropy parameter 7 (we may assume < 7 < ir) is related to the 
quantum group deformation q by q = exp(?7). This R-matrix (2.17) is unitary for 



12 



real 9. Crossed unitarity (2.5) holds for 77 = ry . It is also convenient to work with 
slightly modified local vertices 

[t ab {9)\cd = R b ct(8 ~ <7/2) (2.18) 

in order to construct the row-to- row monodromy matrix T a fr(#, Co) . For this R- matrix 
the diagonal solutions K ± (9) of eqs.(2.12) and (2.15) turn to be[8,9] 

K ± (9) = K(9,Z±) 

where 




K(9,S) = ——[ (2.19) 
smh£ ^ sinh(£-#) J 

Here and are arbitrary numbers that parametrize this boundary conditions 
compatible with integrability. (For the general solution of eqs.(2.12) and (2.15) in the 
six- vertex model see ref. [11].) 

We are interested on boundary conditions yielding a quantum group covariant 
framework. [SU(2) q for the six- vertex model]. For periodic boundary conditions the 
quantum group transformation properties of T a fr(#,cj) and t(9,u) are not simple [15]. 
A SU{2) q invariant XXZ hamiltonian requires fixed boundary conditions and special 
end-point terms[6]. 

The XXZ hamiltonian follows from dt(9,uj = 0)/d9 evaluated at 9 = . One 
finds from eqs.(2.14),(2.17), and (2.19) [8] 

= ~oJ2 + a > V n+i + cos 7 « +1 ) (2-20) 

+ - sinh(i7)(crf coth^_ + a z N coth^ + ). 

As one can check, the quantum group invariant case corresponds to £± = ±00. We 
choose therefore £± = ±00 in order to built a SU(2) q covariant framework for the 
six-vertex model and its scaling limit. This choice corresponds physically to boundary 
conditions of Dirichlet type. In that case 

(<*» \ 

K ± (9)=eMT9cr z )=[ q ^ \ (2.21) 

Now, in order to find the SU(2) q content of the YB algebra (2.14) and its associate 
Bethe Ansatz construction, we start by computing the 9 — > 00 limit of the U a b(9,oj) 
operators. 
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The 6 — > ±00 limit of the row-to-row monodromy matrix 



T(9,u) = 



A{9) B{9) 
C{9) D{9) 



yields SU(2) q generators [1,15]. We have 

*-±°VFiJV7/2 exp(±2 7 J 2 ) [l + 0(e^)] , 
D ^9^±oo eTiNj/2 exp( ^ 7 j z ) [i + ( e ^)] , 

^ e^±oo ± 2eTi JV 7 /2 T e sinh (, 7 ) j_ (t ~ )T7) [ 1+O ( e =F20)] ) 



where 



N k-l . N 

fe=i j=i /=fc+i 



and 



1 N 

Jz = oEw- 



Notice that J + = J+(c<j, 7), J_ = J -{Co, 7) and J 2 are 5C/ (2) g generators with c 
, for they obey the commutation rules 



sin(27 J z ) g 2J * - g 2J * 

i J +' J -J = — — = — : — :=i — 



sin 7 q — q 

[J Z ,J±]=±J± 

For the boundary conditions £± = ±00 , one sets 



U{9,Cj) 



A{9) B{9) 
C{9) V{9) 
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From eqs.(2.10) and (2.21) it follows that 



A{0) =e^' 2 ~ e B(9) C(-9) 

B(6) =e~ i ^ 2+0 A(9) B{-6) 

C{6) =e^ /2 - e D{9) C(-6) 
V{9) =e- l ^ 2+e C{9) B{-6) 

Let us compute the 9 — > oo limit of these operators. We find from eqs.(2.23) -(2.28) 

A(9) 6 =°°- e -«TW exp(2^) [1 + 0{e- 2e )] , 

B(9) d =°° - e-^ /2 sinh(* 7 ) J_ [1 + 0(e- 29 )} , 

C{8) 6 =°° - e-W sinh(* 7 ) J+ e^ J * [1 + 0(e~ 2e )} , (2.29) 

V{9) ° =°° 2e-^' 2 - e sin [(J + J z ) 7 ] sin [(J - J z + 1) 7 ] 

The operator J is defined through the SU(2) q Casimir invariant C q 

r l (i T T t \ i ^ an/ sin2 (7Jz) sin 7 (J + 1) sin 7 J 

C 9 = n( J + J - + J J+) + cos 7 = — (2.30) 

£ sin ( 7 ) sin 7 

As one can see from eqs.(2.29) the asymptotic form of U a i,(X,u) is related with the 
SU{2) q generators in a very clean way. 

Let us now consider the transfer matrix t(\,Cu) [eq.(2.14) ]. We find when £± = 

±oo 

t(8, u) = e -^l 2 - 6 A{9) + e^/ 2+e V(8) (2.31) 
Now, when 9 — > oo , inserting eq.(2.29) in (2.31) yields 

t(0, a)) * =°° 2 cos[ 7 (2 J + 1)] + 0(e- 26 ) (2.32) 

That is, the asymptotics of the transfer matrix expresses up to numerical constants 
in terms of the q-Casimir C q [eq.(2.30) ]. 



e -i~t/2+e A ( Q j D (_ 6 j 

- e^ /2 - e B(9) A(-9) 

(2.28) 

- e^l 2 - 6 D{9) A(-6) 
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The operators U a b(9, Co) fulfil the Yang-Baxter algebra (2.13) with R-matrix (2.17) 
. Letting 9' — > oo in eq.(2.13) we can compute the relevant commutation of U a i,(0, Co) 
with the SU(2) q generators. We find 

[A(9), J_] = -B(0) e^ h+e ' [V(9), J_] = B{9) e i T( J «+ 1 )+ fl (2.33) 

It is well known that 5(#) and C(#) act as lowering and raising operators for J z , 
while A{9) and D{9) commute with it [1] . We find here the same properties for the 
elements of U (9, Co), 

[B{9), J Z ] = B(9) , [C{9), J z ] = -C(9) , 
[ A(9), J z \ = [ V(9), J z \ = [ t(9, Co), J z \ = 

In addition, we find 
e -*7J« B(9), jJ\ = , B{9) J+ = J+ B{9) + [ e ~^ V(9) - e~ e A{9)} 

(2.35) 

Using now eq.(2.31) yields 

[t(6,Q), J-} = (2.36) 

One can analogously prove that 

[t(8,u>), J+] = (2.37) 

Therefore, the transfer matrix t(9,uj) is SU(2) q invariant. As a corollary, we see that 
Hxxz (see eq.(2.20)) is S'6 r (2) a -invariant, as shown in ref. [6] by direct calculation. 

We investigate in the next section the Bethe Ansatz construction of eigenvectors of 
t{9, u). The first consequence of the the SU (2) 9 invariance of t(9, Co) is the degeneracy 
of its eigenvalues with respect to the quantum group. 
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3. The Bethe Ansatz and the SU (2) q group 



In section II we developped the Yang-Baxter framework with boundary conditions 
adapted to the SU (2) q invariance. We call this a quantum group covariant framework. 
Let us now investigate in such scheme the eigenvectors of the transfer matrix t(6,uj) 
in the algebraic Bethe Ansatz [8]. These eigenvectors can be written as 



*(v)=B(v 1 )B(v 2 )....B(v r )n (3.1) 
where Q is the ferromagnetic ground state and v = (v i, V2, ■ ■ ■ , v r ). 



Q= 1 oWo) (o 1 (3 ' 2) 



and 



^inh 9/9 

® {9) = sinh(2fl - Ty) ^ m = ^ m A{ ~ 6) ~ ^ B{ ~ 9) m (3 ' 3) 



Notice that use has been made of the Yang-Baxter algebra (2.8) to obtain the last 
form in eq.(3.3) from the expression (2.28) for B(9). The numbers i>i,i>2, ■■■■,v r , (0 < 
r < N/2) must be all distinct roots of the set of algebraic equations 



-q smh[v m -uj k + «7/ 2 ] sinh[w m + u k + ry/ 2 ] 



k 



L = A sinh[w m -uj k - ry/2] sinh[v m + u k - ry/ 2 ] 

-A- sinh[t; m - v k + ij] sinh[t; m + v k + ij\ (3.4) 
k =it^m sinh \- Vm ~ v k- h] sinh[w m + v k - ij] 
1 < m < r 



These Bethe Ansatz Equations (BAE) guarantee the vanishing of the so-called un- 
wanted terms. That is, vectors in t(9,uj)^(v) which are not proportional to ty(v) . 
They arise, together with the wanted terms, from the repeated commutation of t(6, u) 

with B{vi)B{v 2 ) ■ ■ ■ 

B{v r ) by means of the Yang-Baxter algebra (2.13) . Let us point out that the order of 
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the factors B(vj) in eq.(3.1) is irrelevant since the Yang-Baxter algebra (2.13) implies 



B(9), B{6') 



(3.5) 



From eq.(3.3) it is evident that B{6) is an odd function of 



B( 



-B{6) 



(3.6) 



Consequently, a direct check shows that the BAE (3.4) are invariant under negation 
of any single unknown v j . This implies that it is enough to consider BAE roots 
with strictly positive real parts. In particular, purely imaginary pairs (±iry) are ruled 
out. Moreover, one verifies by explicit calculation that B(i7i/2)Q = . So that the 
selfconjugate root v = i9/2 is also ruled out. 

The eigenvalue of t(9,u) on ^/(v) is given by 



A(9;v) =A+(0;v)+A-(0;v) 



(3.7) 



N 



sinh(26 l ± ry) -i-r sinh[# — Uk + «t/2] sinh[# + Uk + «7/2] 
^ ' V ' ~ sinh(2# + «7) Al s inh[^ - cj fc - 17/2] sinh[# + cu fc - 27/2] 

-pj sinh[6> - v m =F ry] sinh[6> + t; m =F ry] 
J- A sinh[# - v m ] sinh[# + t> m ] 

Taking logarithms, eq.(3.8) become 



(3.8) 



N 



^2 [<t>( v m ~ ^h, 7/2) + (f){v m + CJfc, 7/2)] 



27r/ m + [4>{v m - Vk, 7) + 0(^m + ^fc, 7)] 

1 < m < r 



(3.9) 



where the I m are positive integers[16] and 



sinh(ix + X) , .„ . 

0(A,x =ilog — -4 with 0(0, x) = 0. 

smh(zx — A) 



(3.10) 



Before solving these equations in the N — > 00 limit, let us study the — > 00 behavior 
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of A(9; v) in order to match with the asymptotics (2.29) . We find from eqs.(3.7)-(3.8) 



A(9;v) 



2cos[7(AT + l-2r)] [1 + 0{e~ M )} 



(3.11) 



Here all roots i>i,i>2, ....,v r are assumed to be finite. Actually this must be the case 
unless 7/V is rational (see sec. 5). Comparison of eqs.(2.32) and (3.11) shows that 
there is only one positive (or zero) solution 



except when 7/7r is rational where other possibilities may happen. In addition, we 
have from eqs.(2.34) and (3.2) that 



Eqs.(3.12) and (3.13) show that J = J z for Bethe-Ansatz eigenvectors. That is, they 
are maximal weight vectors for the quantum group for non-rational values of 
7/V . Thus, 



In addition, the kernel of J + is stable under variations of 7 . Therefore eq.(3.14) 
holds even when 7/7? is rational (see below). Each eigenvalue A(9;v) has (at least) a 
degeneracy (2J+1) for t(9,uj) , since the vectors 



are all eigenvectors of t(9,u) with the same eigenvalue thanks to eqs.(2.36)-(2.37). 

Of course, when 'j/n is rational J may not be equal to J z for some BA states 
as follows from eqs.. (2.32), (3.11) and (3.13). 

Let us give an algebraic proof of the highest weight condition (3.14) for the BA 
states. We shall apply J+ to *f?(v) given by eq.(3.1) permuting J + through the 




(3.12) 




(3.13) 



J+ty(v) = 



(3.14) 



J_*(i5), (J-) 2 V(v),..., (J-) 2J *(v), 
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B( v j), (1 < i < r ) with the help of eq.(2.35) . Finally we shall use that 

J+tt = 

Actually it is more convenient to use the operator T>{9) 

1 



(3.15) 



sinh(2# - ij) 
Eq.(2.35) can now be written as 



[V(0) sinh29 -A(9) sinh(i 7 )] 



(3.16) 



J+ B(9) = e t7 B(9) J+ + [e-° A(9) - e e V 
We find from eq.(3.1) after using r times eq.(3.17) : 

T 

J+V(v) =^exp[i 7 (^^ -r + 2j)} B( Vl )B(v 2 ) B( Vj ^) 



(3.17) 



(3.18) 



[e- Vj A(vj) - e^P(uj)] B(y j+1 )....B(y r )Vt 



where we also used eqs.(3.13) and (3.15) . Now we can push the operators A(vj) and 
T>(vj) in eq.(3.18) to the right using the Yang-Baxter algebra (2.13). That is, 

m *w = sinh(g "wr^ sin !!S! + l" 27) m a{6) + 



sinh(# - 9') sinh(fl + 9') 



sinh(2# — i'-f) sint^ry) 
sinh(2#) 



B{9) 



A{9') 



sinh(# - 9') sinh(fl + 9') 



(3.19) 



sinh(# - 6') sinh(# + 9') 



sinh(2# + i'-f) sinh(i7) 
sinh(2#) 



B{9) 



A(9>) 



V(9') 



sinh(fl + 9') sinh(# - 9') 



We find from eqs.(3.18)-(3.19) an expression of the form 



J+V(v) = c j(v) B( Vl )B(v 2 ) B{v j - 1 )B(v j+1 )....B{v r )n (3.20) 

where the Cj(v) are c-number coefficients. The easier is to first compute c\{v) and 
then deduce the rest by symmetry of the arguments (cfr. eq.(3.5) ) . c\{v) follows by 
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permuting in the first term of eq.(3.18) the bracket e Vj A(vj) — e Vj V(vj) through 

B{v2)....B{v r ) using eqs.(3.19) and keeping just the first terms in eqs.(3.19) . That 
is, those proportional to B(vj) A(v\) and to B(vj) T>(v\), (2 < j < r). Collecting all 
factors, it is easy to see that c\{v) and the Cj(v) (1 < j < r) are proportional to the 
BAE (3.4) and hence identically zero. We have therefore proved that all BA vectors 
are highest weights [3,17,18]. 



4. The quantum group covariant ligth-cone approach. 

The light-cone approach is a direct way to give a field theoretical interpretation to 
a lattice model and furthermore obtain its scaling limit as a massive QFT. We start 
from a diagonal lattice that we interpret as a discretized two-dimensional Minkowski 
spacetime. The sites in the light-cone lattice are considered as world events. Each 
site has then microscopic amplitudes associated to it which describe the different 
processes that may take place and must verify general properties like unitarity. 

The light-cone approach with periodic boundary conditions has been investigated 
in refs. [19-20]. Here we consider this approach in the case of fixed boundary conditions 
leading to a quantum group invariant framework. For the six-vertex model, this can 
be done as follows. Consider the SU{2) q invariant transfer matrix t(6,uj) constructed 
in sec. 2, and identify the inhomogeneity parameters u>i, u>2, <^jv with the cutoff 
rapidities of the diagonal lattice. Namely set 

ujj = (-iy e (4.i) 

Then define 

We shall now show that U (O) has the natural interpretation of the unit time evolution 
operator on the diagonal lattice. To this end we introduce some convenient compact 
notation. We denote by the auxiliary horizontal space and label by j, (1 < j < N), 
the vertical spaces. Then 

L J (6)=R 0j (6-uj J )P 0j (4.3) 

is a local vertex acting only on the jth. vertical space. Here Rjk(0) stand for the 
six- vertex R-matrix (2.17) acting on the spaces j and k , (0 < j, k < N) and is 
the usual exchange operator. In practice, the ab matrix element of Lj(8) in the 0th. 
space act on the jth. space as the operators t^iO — cuj) defined in sec. 2. In terms of 
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the Lj{9) we can write the transfer matrix as 

t(e+ij/2,u) = Tr {K+(9 + il )L N (9) Li(9)Kq (9)L 1 (-9)~ 1 L N {-9y 1 } (4.4) 

Since the R-matrix is regular, R(0) = 1 , we have 

Lj(6) \e^ = Poj (4.5) 
Hence, as soon as uoj = (— 1) J 6 , one obtains 

^2fc(@) = PQ,2k , L 2k +i(Q) = S 0j2 k+i (4.6) 
where Sjk = Rjk{2Q)Pjk ■ Inserting eq.(4.6) into eq.(4.4) yields for odd N 

t(e + i-y/2,uj = 0) 
=Tr {K^(9 + i'y) S 0j nPo,n-iSo,n-2Po,n-3 Sq^Pq^Sq^Kq 

Pq,\Sq,2Pq,2, Sq^n-iPq,n} 

=Tr {K£(6 + ij) S 0jN S N -2,n-i Si j2 Po,n-iPq,n-3 Po,2K (<~>) 

Po,iPo,3 Po,nS 2 ,3 SW-i,jv} (4.7) 

=Tr {K$ (9 + ry) S 0jN P 0jN } S u S n _ 2 ,n-iPo,nPo,n-i Po,2 

Po,iPq,s----Po,nK^(Q)S 2 s Sjv-i,jv 

=Tr {K+(9 + t 1 ) Ro,n(2&)} Ru(2G) R N - 2 ,N-i(2&) 

1*7(6) i? 23 (26) R N ^ N (2&) 

Actually, the computation up to now is completely general and would hold for any 
R-matrix and K- matrices. Let us exploit now the explicit form of R(0) and K ± (9) 
for the six-vertex model [eqs.(2.17)-(2.19)]. Specialising to the SU(2) q invariant case, 
we obtain 

/q-(e) = g^O)' 1 , Tr {K+(9 + <7)i2o^(26)} = g N (0) ^G^t) ^ 
where 

g{9) = exp[-9a z ] = I \ (4.9) 

Thus 

U(Q) — R12R34 ■ ■ ■ RN-l-e,N-e9 Nd^ -^23-^45 • • • Rn-2+e,N-1+e (4-10) 

where e = [1 — (— l) N ]/2, gj is the matrix exp[0er z ] acting nontrivially only on the 
two-dimensional vector space attached to the j th link, and Rj^ is the 6V R-matrix 
R(2Q) acting nontrivially only on the tensor product of the j th and k th vector spaces. 
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U(0) is clearly the unit time evolution operator on a diagonal lattice with "re- 
flecting" type boundary conditions, g(9) [g^)^ 1 ] acting for collisions on the right 
[left] wall. Graphically U(Q) can be depicted as in fig.6. By taking powers of U(Q) 
one then obtains the evolution on the diagonal lattice for any discrete time (Fig. 7). 
By construction, this time evolution is SU(2) q invariant. 

The eigenvectors of U(Q) can now be written as in sec. 3, namely 



=B(v 1 )B(v 2 )....B(v r )Q 



(4.11) 



where uj = (— 1) J 9 and the Bethe Ansatz equations (3.4) take the form 



i N 



sinh[w m - 6 + ry/2] sinh[w m + 6 + ry/2] 
sinh[f m — — ry/2] sinh[t> m + 6 — ry/2] 

sinh[w m - Vk + i"f] sinh[t> m + Vk + %l\ 



n 



(4.12) 



1 < m < r 



sinh[w m -v k - ry] sinh[w m + v k - 17] 



The associated eigenvalue of [/(©) is given by 

7y / n \ TT sinh [ e - - sinh[9 + t; m - yy/g] 

-"■j 1 - smh[6 -w m + «7/2J smh[6 + v m + «7/2J 

Let us now discuss the problem of unitarity for the time evolution defined by U (0) . 
The R-matrix (2.17) is unitary for real 6. Hence the evolution operator U is unitary 
too, for real 0, up to the boundary effects due to the term g^g^ 1 in eq.(4.10). As 
is evident from figs. (6-7) [eqs.(4.10)] there exist unitarity violating boundary effects 
(factors g(6)N and giO)^ 1 ). Then, the reflection of the bare particles on the bound- 
aries is affected by " leaking" . However, on average this leaking compesates since we 
have opposite factors g(6)N and giO)^ 1 on right and left sides respectively. Thus 
one can expect U(Q) to still have unimodular eigenvalues. Indeed, this is the case 
from eq.(4.13) provided, as it is usually the case, the BA roots are either real, self- 
conjugate (Imv = 7r/2) or organised in complex conjugated pairs. As seen below, 
the exact diagonalization of U(Q), shows that all its eigenvalues are unimodular, so 
that there exists a similarity transformation mapping U(Q) to an explicitly unitary 
operator. On the other hand, it is natural to expect that in the thermodynamic limit 
(N — > 00) different boundary conditions become equivalent (the model has a finite 
mass gap). This suggest that the above mentioned similarity transformation reduces 
to the identity as iV — > 00. We conclude therefore that the light-cone 6V model 
with fixed b.c. described by the evolution operator U of eq.(4.10) is another good, 
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integrability-preserving regularization of the SG model. Its advantage over the more 
conventional setup with periodic b.c. is that U(Q) is explicitly SU(2) q invariant even 
for finite N . 

According to the general light-cone construction the physical lattice hamiltonian 
can be defined in terms of U(Q) as 

if = - log 17(6) (4.14) 
a 

where a is the lattice spacing. By a judicious choice of the logarithmic branches the 
energy eigenvalues can be written from eqs.(4.13)-(4.14) 

E = a- l J2^(vj) (4.15) 

3=1 

where 

„ , /cosh 26 cos 7 — cosh2tA 

e (v) = -7T + 2 arctan ' , 4.16 

\ sinh 29 sin 7 / 

Notice that eo(v) is smooth and negative definite for real v. The specific choice of 
the logarithmic branch in passing from eq.(4.13) to eqs.(4.15), (4.16) is dictated by 
the requirement that eo(u) should correspond to the negative energy branch of the 
spectrum of a single spin wave over the ferromagnetic state Q. 

The physical ground state and the particle-like excitations are obtained by filling 
the interacting sea of negative energy states. This sea is described by a set of real 
v— roots of the BAE with no "holes" (in the Appendix we give a detailed exposition 
on the treatment of the BAE, also to clarify some rather subtle matters). Excitations 
correspond to solutions of the BAE which necessarily contain holes and possibly 
complex v— roots. The crucial point is that in the limit N — > 00 only the number 
of real roots of the sea grows like N, while the number of holes and complex roots 
stays finite to guarantee a finite energy above the ground state energy. Hence these 
solutions of the BAE can be described by densities p(v) of real roots plus a finite 
number of parameter associated to the positions x\,X2, ■ ■ - x v of the v holes and to 
the location of the complex v— roots. 

The continuous theory is defined by the double limit a — > 0, 9 — > 00 taken in 
such a way that a finite mass gap m emerges. The explicit result in the periodic case 
which will be shown to hold also here is [1,19,20] 



m = lim 

a^O 6^00 



— e a 
a 



(4.17) 



[The double limit is taken such that the l.h.s. is finite and non-zero]. The invariance 
of £7(9) under the quantum group SU(2) q allows to perform a transformation of the 
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original vertex model where states are assigned to the links and interactions as 'bare' 
scattering amplitudes, to the vertices, into a diagonal lattice face model, where states 
are defined on the plaquettes. The plaquettes are the sites of the dual lattice and the 
interaction involves the four plaquettes around each vertex, or equivalentely, the four 
sites around each face of the dual lattice. 

On the faces of the diagonal lattice one introduces positive integer-valued local 
height variables £ n (t), t according to (see fig. 2) 

e (t) = 1; hit) = 2 

W ' W (4.18) 
ln+i(t)=ln(t)±l,n = l,2,...N-\ 

The configurations {h(t), ■ ■ ■ ^at(O) a ^ an Y fixed discrete time t are in one-to-one 
correspondence with the SU(2) q multiplets of the 6V model of dimension £ n (t) = 
2 J + 1 (notice that £jv(t) can be chosen to be time-independent thanks to SU(2) q 
invariance). Now the matrix elements of U between these highest weight states define 
the unit-time evolution in the face language. 

This vertex-IRF correspondence is well known : in practice, it amounts to an 
application of the q-analog of the Wigner-Eckart theorem. Indeed, by SU(2) q invari- 
ance, the matrix elements of U(Q) in the subspace with definite total q-spin J = j , 
for 7/7r not a rational, can be written 

<jm;r\U(e)\j'm';r' >= 5 jf S mm > < r\\U(Q)\\r' >j (4.19) 

where m = —j, —j + 1, +j labels the J z eigenvalues and r the degeneracy of the 
J = j subspace which has dimension 



Vf ] (4.20) 



< r||f/(@)||r / >j thus denotes the reduced matrix elements, and a specific choice of 
basis is required to give explicit expressions for them. The most natural basis, which is 
the basis useful for the vertex-IRF correspondence, is that obtained by the successive 
composition of the q-spin 1/2 basic constituents assigned to the links. Consider a 'time 
zero' line cutting the diagonal lattice as in fig. 8. By intersecting N links, it identifies 
the Hilbert space of the vertex model as (®[l/2])^ , where [j] denotes an irreducible 
representation of weight j . On the other hand, the line passes also through a well 
defined set of plaquettes, including the half-plaquette at the extreme left and right. 
Then, a set of configurations of a local height variable I , assigned to each plaquette, 
can be constructed as follows. Assign I — Iq — 1 to the half plaquette on the left, the 
0th. plaquette. Then pass to the next one on the right , the first plaquette, following 
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the time zero line. This cuts a link carrying a spin 1/2 representation. So that, at this 
stage the Hilbert space is just the representation [1/2] itself, i.e. J = 1/2 . We then 
set l\ — 2 J + 1 = 2 . Following the line, we now cut another spin 1/2 link and arrive 
at the second plaquette. This Hilbert space is now the direct sum [J = 0] © [J = 1] 
and the local height 1 can take two values, 1% = 2j + 1 = 1 or 3 . By repeating this 
procedure n times (with n < N), we land in the nth. plaquette after having cut n 
spin 1/2 links : the Hilbert space is (<g> [1/2] ) n and contains irreps with J running 
from or 1/2 up to n/2 . Hence the local variables In take the values 



By construction, \ln — ln-\\ = 1 , since the In follow the composition laws of SU(2) q 
representations 



which for q = exp(ry) not a root of unity, are just those of the usual SU(2). 

When n = N we arrive at the half-plaquette on the right. Here the values of In 
are just the dimensions of the irreducible representations in which the full Hilbert 
space (<g)[l/2]) iV can be decomposed. All together, a specific choice of fo, h, - ■ ■ ,In ( 
Iq and l\ are fixed to Iq — 1, l\ — 2 by construction) defines one and only one of the 
Fj(N) irreps with J = j. , if In — 2j + 1 . In other words, we have constructed a 
map from the N + 1 plaquettes at 'time zero' and the Bratteli diagrams giving the 
composition rules for the tensor product of N spin 1/2 reprersentations (see fig.9). 
We are now in position to identify the degeneracy label r with a specific path in the 
Bratelli diagram, i.e. r = (Iq, h, ■ ■ ■ ,In) ? with Iq — 1 and In — 2j + 1 . 

It remains to evaluate the reduced matrix elements 



where the height variables l'o, . . . , are associated to the plaquettes crossed by 
the 'time one' line (see fig. 8). Notice that Iq — l'o — 1 by definition, while In = 
I'n — 2j + l by SU(2) q invariance. However, in spite of the factorised form of U(0) 
[eft. eqs.(4.10)], these matrix elements cannot be written in a factorised form, with 
each factor depending locally on the variables Iq, h, . . . In and I'q, . . . ,In ■ This 
is so because the full unit time evolution operator U(Q) is SU{2) q invariant, but 
each single i?/ Cj fc + i(20) which enters is not. Only the peculiar boundary terms g\ 
and qn enforce SU(2) q invariance. Nevertheless, it is easy to pass to a different 
vertex representation where the SU{2) q symmetry holds in a local sense. We define 




[J]® [1/2] = [J +1/2] ©[J -1/2] 



< rii^e)!!/ >j = < k, h, ■ ■ ■ , i N \u(e)\i' , i\, 



■ , In > 



(4.21) 



2G 



a similarity transformation on U (0) 

U(Q) = GU{Q)G- 1 (4.22) 

where 

n 1/2 -1/2 1/2 e-1/2 OQ ^ 

G = #i £ 2 ...^jy (4.23) 

l— (— i) N 

e = — 1 2 ' and is exp[— Bcr z ] acting on the j-th. spin 1/2 space. Now comparing 

eqs.(4.10) and (4.22)-(4.23) it is easy to verify that t/(0) can be written, e. g. for 
even N 

U (O) = -R12-R34 • • • -RjV-l,iV-R23-R45 • • • Rn-2,n-i (4.24) 

where 

R k ,k+i = 9 k 1/2 gtZRk,k + im9l /2 g k l{ 2 (4-25) 

One can explicitly verify that Rt,k+i is SU(2) q invariant; that is 



[-Rfc,fc+1, (<Tz)m + (Vz)m+l] =0 
[R kM1 , (<T^ /2 ) m (a±) m+1 + (a ± ) m (q^ 2 )} =0 



(4.26) 



The reduced matrix elements of U (0) can now be expressed in a factorised form , 

< h, h,..., In\u{Q)\i' , i\, ...,i N >= w[wiwi . . . W' N _ 1+e W 2 W 4 . . . W N - e , 

where 

W m = W(l m -l,lm+l\lm,l'm;Q) , W m — W(l' m -i, l' m +l\lm, I'm',®) (4.27) 



and 
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W(x,y\u,v; 0) = 8 UV - ■ ,T n 7 ■ S V ^ ^ <V (4.28) 

smh(2B + z7j [x\ q 

with the standard notation [x] q = (q x — q~ x )/(q — q^ 1 ) = sin^x)/ sin7. This com- 
pletes the transformation from vertex to faces description. The global time evolution 
in the faces (or heights) language is obtained by taking matrix products of 

< b,h,---, In\U(Q)\1'o, I' i,..., In > 

The boundary conditions are Iq = constant = 1 on the left, while In = constant = 

2j + 1 on the right, if the reduction is performed onto the dimensional space of 

irreps of q-spin J = j . With the natural constraint that \l — l'\ = 1 , when I and F 
sit on neighboring faces, the weights given by eqs.(4.27)-(4.28) define the ABF-SOS 
model in the trigonometric regime [13]. 
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5. Analysis of the Bethe Ansatz equations 



We investigate in this section the solution of the BAE (4.12) associated to the 
quantum group covariant BA in the light-cone approach. It is convenient to relate 
them with the BAE for periodic boundary conditions (see for example [1]). Define 2r 
variables Xj as 

Xj = Vj , Xj +r = -v r -j + i , 1 < j < r (5.1) 



Then eqs.(4.12) can be rewritten as 



sinh[A m — 6 + ry/2] sinh[A m + 6 + i r y/2]~\ N sinh(2A m + 27) 

sinh[A m — — ry/2] sinh[A m + — 27/2] J sinh(2A m — ry) 

2r ' (5-2) 

nsmh[A m - A fc + 17] 
— — — , 1 < m < 2r. 
fc=i smh[A m - X k - vy\ 

These equations are like the BAE for periodic boundary conditions on a 2N sites line 
and with an additional source factor 

sinh(2A m + ij) _ sinh(A m + ry/2) sinh(A m - i^g) 
sinh(2A m -i7) sinh(A m - ry/2) sinh(A m + i^p-) 

More important, we have the following constraints on the roots X m : 

a) the total number of roots is even (2r) and they are symmetrically distributed 
with respect to the origin according to eq.(5.1). 

b) Xm = and purely imaginary X m are excluded as roots. 

Let us start by considering the antiferroelectric ground state. It is formed by 
roots with fixed imaginary part equal to n/2 . Therefore, we shift X m — > A m + in/ 2 
and take X m real for the ground state. Now X m = is excluded whereas it may 
be present for periodic boundary conditions (PBC). Therefore, the ground state for 
eq.(4.12) is a one-hole solution of eq.(5.2) with the hole at A m = O.We also redefine 7 
into 7r — 7 in order to agree with the conventions of refs. [1,19,20]. With this choice, 
7 is related to the sine-Gordon coupling constant f3 by (3 2 /(8n) = 1 — 7/71. 

Let us now consider the iV — > 00 limit where a continuous density of BAE real 
roots can be introduced 

" (A '») = i?„ N(X, n+ \ - X m ) (5 ' 4) 

This function must always be even in our fixed boundary conditions case. Taking 
logarithms in both sides of eq.(5.2) and using eqs.(5.3) and (5.4) the BAE yield linear 
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integral equations through the usual procedure 



N [$'(A + 0, 7 /2) + <t>'(A- 0,7/2)] + $'(A, 7/2) - $'(A, ^-Jl) 



+00 

!,, n(,A ch'(\ — _i_ 

(5.5) 



= 27iNp(\) +N J dnp{n)&{\-n,i) + 



2N C 2M h 

[S'(A - z u 7) + $'(A - ^,7)] + 27nf (A) + 2vr ]T 5(A - 9 h ) 
1=1 h=l 

where we assume 2N C pairs of complex roots (z^, Zk) and 2M^ holes 6^ , symetrically 
distributed with respect to the origin. The roots density p(A) is connected with cr(A) 
, the derivative of the counting function through 



1 , 2M h 

P(A) = a(A) - -5(A) - - £ 5(A - 6l ft ) (5.6) 



The Fourier transform solution of eq.(5.5) gives for the ground state 



+00 

<to(A) = ! ^ e* fcA flro(fc) 



(5.7) 



_ ... cos(fce) + l/(2iV) 1 sinh[A;(7r-37)/4] 

= 7-7^ 77TT h 



cosh(A;7/2) 2N cosh(k-f/2) sinh[A;(7r - 7) /4] 



Notice that the first term is twice the PBC ground state density (eft. ref.fl]). The 
second term is the boundary effect and leads to a A^ 1 surface correction to the free 
energy [12]. 

The densities of real roots for excited states, that is in the presence of (non-zero) 
holes and complex roots are analogous to the PBC case. We find from eq.(5.5) : 



(7(A) = (70(A) + l[<7 h (A) + <7 c (A)] (5.1 



Here <7/j(A) and er c (A) stand for the holes and complex roots contributions, respec- 
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tively. The explicit expressions of these densities are: 



2M h 2N C 

<7 h (A) = J>(A - h ) , a c (\) = J2Q^-^Vl) (5-9) 
h=l h=l 



where Zk = o~k + i"Hk with % > and 



+00 



_ / dfc sinh(7r/2- 7 )A; ifcA 

AJ 7 2tt 2sinh(vr-7)A;/2 cosh( 7 A;/2) 

— oo 

Q(X, rj) — - p(X - irj) - p(X + irj) , r] < 7 < tt/2. 



(5.10) 



Here 



^(A, 77) = $ ( ) ( 5 - n ) 



7T 7T 



Complex roots with Im (z) > 7 are called wide pairs whereas roots in the region 
Im (z) < 7 appear as two-strings (Imz = ±7/2) or in quartets[14]. 

Let us now derive the higher-level BAE. We follow a procedure analogous to 
ref.[14] for the PBC case. The only novelty with respect to the PBC case is the extra 
term — $'(A, (ir — 7)/2) in the source and the obligatory presence of a hole at A = . 
In order to derive the higher-level BAE one has to compute 

cr (A) + <7 (A - ry) 

We find from eq.(5.7) 

<7o(A) + a (X - 27) = ^Jx^( 2X - *7, 7) (5.12) 

The r.h.s. comes entirely from the second term in eq.(5.7) . Notice that the l.h.s. of 
eq.(5.12) exactly vanishes in the PBC case. 

We define the parameters Xj as m the PBC case [14] : 

X = real part of the two-string position, 

X = z c - ry/2, z c + ry/2 for a quartet (z c , z c , z c - 27, z c + 27), 
X = z w - 27/2, z w + ry/2 for a wide pair (z w , z w ). 



30 



Then, the BAE take the following form in the N = oo limit 



tt sinha[xj - 9 h + h/2] _ 

M P (5-13) 
sinha[xj-X*-»7j 

where a = 1/(1 — 7/77) and the use of eq.(5.12) yields the phase factor e l ^ 2xj ^ 
in the r.h.s. . The v = 2M/, holes and the M p complex pairs are all distributed 
symmetrically with respect to the Re A = axis. Hence, restricting to positive real 
parts, we have 

tt sinh a[xj - Oh + jr/2] sinha[xj + flfr + yy/2] _ 
Al sinh a[xj -9 h - ij/2] sinh a[xj + 6 h - ij/2] 

M yi sinh a[xj - Xk + jr] sinh + Xk + jr] ( 5 ' 14 ' ) 
fe= '|^ sinh - Xk ~ *7] sinh + Xfc - *7] 

1 < J < Mp/2. 

To conclude, let us come back to the BAE (4.12). They can be rewritten in a mani- 
festly algebraic form by introducing the variables zj = exp(—2vj) , (1 < j < r). 

zjw -q zj - wq \ N _ yr Zj - z k q 2 zjz k - q 2 
Zjwq-lzjq-wj Al z . q 2 _ Zk ZjZk q 2 - 1 

Mi 

where u> = exp(— 20) and g = exp(ry) is the quantum group deformation parameter. 
According to the general discussion following eq.(3.6), we can restrict the search for 
solutions to eq.(5.15) strictly within the unit circle \zj\ < 1 . Therefore, to any 
unordered set of distinct numbers fulfilling eq.(5.15), z\, . . . , z r , with \zj\ < 1 , there 
corresponds one BA eigenstate of the transfer matrix (4.4) and hence of the evolution 
operator on the light cone-lattice. Special attention should be paid to the possibility 
that one or more roots zj lay exactly at the origin. This corresponds to Re Vj — > 00 , 
and therefore to the reduction of the corresponding B(vj) to a multiple of the lowering 
operator J_ of SU(2) q [eft. eq.(2.29)]. 

B(oo) = (1 -q 2 )J- (5.16) 

In effect, the point at infinity constitutes the only exception to the requirement that 
the roots v\, . . . , v r be all distinct. However, this is a very special possibility. Indeed, 
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suppose that a given root in eq.(5.15) , which we can always identify with z\ , lays 
at the origin z\ — . We immediatly obtain from eq.(5.15) : 



q ™ = g 4(r-l) ^ g 4(J z + l) = 1 (51?) 

that is, q must be a root of unit and 7 a rational multiple of n . This correspond 
to the special cases when certain irreducible representations of SU(2) q mix into type 
I reducible but indecomposable representations. This interesting phenomena is dis- 
cussed in sees. 7-9 within the Bethe Ansatz framework in connection with the RSOS 
models. 

Eqs.(5.14) can also be written in algebraic form: 

TT jj w h ~ g jj - W h q = yr Q ~ £fcg 2 ~ f ^ ^ 

f =1 ij w M - 1 iA - w h U ijq 2 - & ijikq 2 - 1 



where 



e - 2 «» W/i = e -2a^ ; $ = e i7 = e <a7 a = (5.19) 

7T — 7 



These new BAE involve only a finite number of parameters and are formally identical 
to the "bare" BAE (4.12)-(5.15), with the holes acting as sources in the place of the 
alternating rapidities ±0. Moreover, from (5.18) one reads out the renormalization 
of the quantum group deformation parameter 

7T7 

q — * q i.e. 7—^7= (5.20) 

7T — 7 

As we shall see later on, this renormalization has a nice physical interpretation for the 
SOS and RSOS models related to the 6V model by the vertex-face correspondence. 
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6. Higher-level Bethe Ansatz and S-matrix 



The holes in the sea of real BAE roots are the particles of the light-cone 6V 
model. They are SU(2) q doublets and can be present in even (odd) number for N 
even (odd). Consider first the even N sector, setting N = 2N' . The energy of a BA 
state with an even number v of holes located at x\,X2, ■ ■ ■ x v can be calculated to be, 
in the N' — > oo limit, 

E = EQ + a~ 1 Ye(x h )+0(a- 1 N~ l ) , e(x) = 2arctan ( cosh W^ ( 61 ) 
^— ^ \sinh7rB/7/ 

where Eq is the energy of the ground state, that is the state with no holes. Eq is 
of order N and is explicitly given in the appendix along with some detail on the 
derivation of eq.(6.1). Now suppose N odd, with N = 2N' — 1. To compare this 
situation with the previous one, we need to slightly dilate the lattice spacing a to 
a' = 2N'a/(2N' — 1), in order to keep constant the physical size L = Na of the 
system. Then eq.(6.1) remains perfectly valid, as A^' — > oo, also for the case of N 
odd, with the same ground state energy Eq. The only difference is that now v is odd. 
Hence, alltogether, we obtain that the number of holes can be arbitrary (unlike in the 
treatment with periodic b.c.) and that the total energy, relative to the ground state 
and in the L — > oo limit, is the uncorrelated sum of the energy of each single hole, 
independently of the complex pair structure of the corresponding BAE solution. The 
BA state with one hole has J = N/2 - M = (2N' - l)/2 - (N' - 1) = 1/2 and is 
therefore a SU (2) q doublet. The 2 V polarizations of a state with v holes are obtained 

by considering all solutions of the higher-level BAE (5.18) with < M < u/2. We 
see in this way that the holes can be consistently interpreted as particles. 

As a lattice system the light-cone 6V model is not critical. The (dimensionless) 
mass gap is the minimum value of the positive definite e(x), the energy of a single 
hole, that is 

e(0) = 2 arctan ( — — '— — ) (6.2) 
\sinh7r9/7 / 

This gap vanishes in the limit — > oo. Hence the continuum limit a — > can be 
reached provided at the same time 6 — > oo in such a way that the physical mass 

m = a _:L e(0) ~ 4a~ 1 e~ 7r9/7 (6.3) 

stays constant [19]. In the same limit we obtain the relativistic expression 

a~ 1 e(x) — > m cosh irx/'y (6.4) 

so that irx/'-f is naturally interpreted as the physical rapidity 9 of the hole. This is 
consistent also in our fixed b.c. framework, provided the limit L — > oo is taken before 
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the continuum limit. Indeed the total momentum becomes a conserved quantity in the 
infinite volume limit and its eigenvalues can be expressed in terms of the BA v— roots 
exactly as in the periodic b.c. formulation. Then one finds that the momentum of a 
single particle takes the required form msinh# in the continuum limit. 

The above analysis shows that a relativistic particle spectrum appears in the 
a — > limit above the antiferro magnetic ground state. For 7 > n/2 one can show 
also that bound states appear, associated to appropriate strings of complex roots, just 
as in the periodic b.c setup. In the sequel we shall anyway restrict our analysis to 
the repulsive 7 < it/ 2 region, where the only particles are the holes, to be identified 
with the solitons of the SG model. Of course at this moment, since the infinite 
volume limit is already implicit, the particles are in their asymptotic, free states: the 
rapidities 6h = n^h/l ma J assume arbitrary continuus values and the total excitation 
energy is the sum of each particle energy and does not depend on the internal state 
of the particles. The situation changes if we consider L very large but finite, since in 
this case the hole parameters still quantized through the "bare" BAE 

(4.12). Indeed, by definition the holes are real distinct numbers satisfying 

Z N (x h ;v 1 ,V2,...,v M ) = ; h = l,...,v (6.5) 

where Z^[x ; v\, . . . , vm) is the "counting function" defined in the Appendix and the 
positive integers I\, . . . , l v are all distinct from the integers ii, . . . , I r labelling the 
r (r < M) real v— roots of the BAE. For N very large, J = N/2 — M finite and 
x < (7/V) In N, the counting function can be approximated as 

Z N (x ; vi, . . . , v M ) = Zoo(x) + N~ 1 F(x ; x\, . . . , x M \ Xi, ■ ■ ■ , Xm) + { N ~ 2 ) ( 6 - 6 ) 

where (x) is the ground state counting function at the thermodynamic limit 

/ s / sinh7rx/7 \ 

Zoo X = 2 arctan — - '/ 6.7 
\ v cosh7rB/7/ 



and 



F(x ;xi,..., x M ;Xi, Xm) = -» lo g Yl S o(n( x ~ x h))So{^(x + x h )) 

M ( 6 - 8 ) 

yr smh g(x - Xj + nyg) smh a(x + Xj + »t/2) 
pi sinho;(a; — Xj — «7/2) sinho;(a; + Xj ~ «7/2) 

In the last expression, the numbers Xj are the roots of the higher level BAE (5.18), 
(5.19), while Sq(9) coincides with the soliton-soliton scattering amplitude of the SG 



34 



model 



Sq(9) = exp i 



dk sinh(7r/27 — l)k sin k9/n 
k sinh(7r/c/27) coshA;/2 



(6.9) 



under the standard identification 7/ it = 1 — /3 2 /8n. 

Combining eqs. (6.5) and (6.6), and taking the continuum limit one obtains the 
"higher level" expression 



M 



exp (-imLsmh9 h ) = J] S (9 h - 9 n ) S (9 h + 9 n ) J\ 



n=l 



l Ji tj w hQ - 1 tj- w h q 



(6.10) 



where, according to (5.19), £j = e~ 2aXj and = e~ 2 ^® h / 7T . Together with the higher 
level BAE (5.18), this last equation provides the exact Bethe ansatz diagonalization of 
the commuting family formed by the v renormalized one-soliton evolution operators 
(see fig. 3) 

Uh =S h , h -i(9 h - . . . S hA (9 h - 9i) g h {-9 h ) S hA (9 h + 9^... S hjl/ (9 h + 9 v )x 
9h{9 h ) S hjV (9 h - 0„) . . . S hjh+1 (9 h - 9 h+ i) 

(6.11) 

where gh{9) = exp 6a and S(9) is the complete 4 x 4 SG soliton S— matrix in the 
repulsive regime 7 < 7r/2 (for brevity we set 9 = ^9 /{-n — 7)): 



S{9) = S (9) 



(\ 0\ 

b c 

a i 

\0 1/ 



S = 6(fl,7) 



c = c(0,7) 



sinh # 



sinh(i7 — #) 
sinh 27 

sinh(i7 — 9) 



(6.12) 



The operators are the values at 9 = 9h of the fully inhomogeneous Sklyanin-type 
transfer matrix T(9 ; Q\, . . . , 9 V ) constructed with S(9 — 9^) as local vertices. Let 
us stress that the higher-level Bethe Ansatz structure just described follows directly, 
after specification of the BAE solution corresponding to the ground state and without 
any other assumptions, from the "bare" BA structure of the light-cone 6V model. In 
particular, this provides a derivation "from first principles" of the SG S— matrix. We 
want to remark that S(9) is also the exact ^—matrix for the elementary excitations 
of the 6V model on the infinite lattice: it is a bona-fide lattice S— matrix. Of course, 
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in this case 9 = Q\ — 62 is the difference of lattice rapidities, which are related to the 
energy and momentum of the scattering particles through the lattice uniformization 



n _i _ . n . . _i / smh^ 

Pj —a ZoohOj/Tc) = 2a arctan — -7— 

\cosh7T0/7 

n _i . . . . ^ _-■ / cosh# 7 - \ 

£/ 7 - =a eHVj/ir) = 2a arctan — - — — J — 
J J \sinh7r9/7/ 



(6.13) 



implying the lattice dispersion relation 

cos aEj 1 2 = tanh(7r9 /j) cos aPj /2 (6.14) 

Thus, we see that the continuum limit only changes these energy-momentum relations 
to the standard relativistic form, (Ej, Pj) = m(cosh#j, sinh^j), without affecting the 
S— matrix as a function of 6. 



7. SOS and RSOS reductions and kink interpretation 

Under the vertex-face correspondence previously described, the light-cone 6V 
model is mapped into the SOS model. At the level of the Hilbert space, this cor- 
responds to the restriction to the highest weight states of SU(2) q : each (2J + 
1)— dimensional multiplet of spin J (for generic, non-rational values of 7/77) is re- 
garded as a single state of the SOS model. In particular, in such a state the local 
height variables £ n have well defined boundary values 1$ = 1 and £n = 2J + 1. 
The ground state of the 6V model, which is a SU(2) q singlet, is also the ground 
state of the SOS model: it has £0 = = 1 and is "dominated", in the thermodi- 
namic limit iV — > 00 by the see-saw configuration depicted in fig. 1. Our boundary 
conditions allow for only one such "ground state dominating" configuration, with 
£n — (3 — (— ) n )/2, while periodic b.c. on the SOS variables i n would allow any 
possibility: l n = {21 + 1 ± (— ) n )/2, with i any positive integer. 

Consider now a BA state with one hole. Then iV is necessarily odd while J = 1/2 
and £n = 2J + 1 = 2. From the SOS point of view this state is "dominated" 
by height configurations of the type depicted in fig.4. Upon "renormalization" , we 
can replace the ground state and one hole configurations with the smoothed ones of 
fig. 5. The hole corresponds, in configuration space, to a kink of the SOS model 
which interpolates between two neighboring vacuum states: the vacuum on the left 
of the kink has i n = (3 — (— ) n )/2, corresponding to a constant "renormalized" height 
£ n = 1, while the vacuum on the right has £ n = (5 + (— ) n )/2, corresponding to 
£ n = 2. Of course, many kink configurations like that depicted in fig. 4 are to be 
combined into a standing wave (a plane wave in the infinite volume limit which turns 
the one hole BA state into an eigenstate of momentum). Let us observe that, by 
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expressions like "dominating configuration", we do not mean that, e.g. the ground 
state, becomes an eigenstate of £ n as iV — > oo. We expect local height fluctuations 
to be present even in the thermodynamic limit or, in other words, that the ground 
state remains a superpositions of different height configurations. The identification of 
a dominating configuration is made possible by the integrability of the model, which 
guarantees the existence of the higher-level BA. In turns, the higher-level BA allows 
us to consistently interpret the ground state or the one hole state as in figs. 1, 4 and 5, 
since the presence of fluctuations only renormalizes in a trivial way the scattering of 
physical excitations relative to the bare, or microscopic, R— matrix (2.17), as evident 
from eq.(6.12). Therefore, we can reinterpret at the renormalized level the vertex- 
face correspondence: the holes, that is the solitons of the SG model, are SU(2) q 
doublets acting as SOS kinks that increase or decrease by 1 the renormalized heights 
£ n . The higher-level BA makes sure that the total number of internal states of v 
kinks interpolating between £ = 1 and £ = 2J is just the number of highest weight 
states of spin J in the tensor product of v doublets 

d " w = fc) " (m - 1) ; m = W\-J P.D 

This parallels exactly the original BA, which provides the 

^=(m)"(m-i) ; M=im~J (7.2) 

highest wight states of N doublets with total spin J. 

Comparing eqs. (2.17), (4.28) and (6.12), we can directly write down the S— matrix] 
of the SOS kinks: 

S(O)*% = S o (O)^6 r . + 6 ij b0,<y) 

It defines the two-body scattering as follows. The first ingoing kink interpolates 
between the local vacuum with £ = i and that with £ = r (r = % ± 1), while the 
second interpolates between £ = r and £ — j (j — r ± 1). The outgoing kinks are 
interpolating between £ — i and £ = s (s — i ± 1, i.e, s = r, r ± 2) and between £ = s 
and £ — j (j = s ± 1). 

When 7 = n/p for p — 3, 4, . . ., the SOS models can be restricted to the RSOS(p) 
models, by imposing £ n < p. In the vertex language of the 6V model this corresponds 
to the Hilbert space reduction to the subspace formed by the so-called type II repre- 
sentations of SU (2) q with q p = —1 [5,6]. In the next section we shall describe in detail 



sm7r sin 7s 
sin sin 77 



(7.3) 
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how the restriction takes place in our BA framework. Here we simply observe that 
the local height restriction, when combined with the finite renormalization 7^7 
of eq.(5.20), provides a strong support for the kink interpretation presented above. 
Indeed, if 7 = ir/p, then 7 = n/(p — 1) and the renormalized heights £ n can take 
the values 1,2, . . .,p — 2. This appears now obvious, since each constant configu- 
ration of £ n corresponds to an £ n — configuration oscillating two neighboring values. 
Moreover under the standard identification of the critical RSOS(p) models with the 
minimal CFT series M p , we see that each light-cone RSOS(p) model has kink excita- 
tions whose S— matrix (7.3) is proportional to the (complex) microscopic Boltzmann 
wights of the RSOS(p — 1) model, under the replacement of 20 = — (—0)= rapidity 
difference of light-cone right and left movers, with (p— 1) times 6 = 61 — 62 =rapidity 
difference of physical particles. This is exactly the pattern found by bootstrap tech- 
niques [21] for the minimal model M p perturbed by the primary operator 0i.3 (with 
negative coupling). 

8. BA roots when q is a root of unity and Quantum Group reduction. 

As previously explained to each M roots solution of the BAE there corresponds 
a highest weight state of the quantum group SU(2) q with spin J = N/2 — M. It 
is well known that when q is a root of unity, say q p = ±1, then {J+) p = {J-) p = 
and the representations of SU(2) q divide into two very different types. Type I 
representations are reducible and generally indecomposable. They can be described 
as pairwise mixings of standard irreps (that is the irreducible representations for q 
not a root of unity) with spin J and J' such that 

\J-J'\<p, J + f =p- 1 (modp) (8.1) 

Notice that the sum of q— dimensions for this pair of reps vanishes. 

Type II representations are all the others. They are still fully irreducible and 
structured just like the usual 577(2) irreps. Since (J±) p = 0, type II representations 
have necessarily dimension smaller than p, that is 

J < ^ (8-2) 

We shall now show how the BAE reflect these peculiar properties of the SU(2) q 
representations for q a root of unity. First of all let us stress that when q is not a 
root of unity (i.e.^/n is irrational), then the BAE cannot possess v— roots at (real) 
infinity. Indeed, since B(oo) is proportional to J_ (eq.(5.16)), a root at infinity means 
that the corresponding BA state is obtained by the action of the lowering operator 
J_ on some other state with higher spin projection J z . But for q not a root of unity, 
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the BA states have J = J z and therefore cannot be obtained by applying J_ on any 
other state. On the other hand, if we assume that one z— root, say z\, of the BAE 
(4.12) lay at the origin, i.e.Revi = +00, then eqs. (4.12) for j — 1 imply 

g 4 ( J+1 ) = 1 (8.3) 

That is, q must be a root of unity. It is now crucial to observe that the remaining 
equations for the non-zero roots (those labelled by j = 2, 3, . . . , M) are precisely the 
BAE for M — l unknowns. This invariance property holds only for the quantum group 
covariant, fixed boundary conditions BAE. It does not hold for the p.b.c. BAE where 
v— roots at infinity twist the remaining equations for finite roots. This twisting reflect 
the fact that the corresponding periodic row-to-row tranfer matrix is not quantum 
group invariant but gets twisted under SU(2) q transformations [1]. Thus, when q is 
a root of unity, BA states (4.11) with one v— root at infinity take the form 

= 00, v 2 , • • .,vm) = (1 - q 2 )J-^(v2,v 3 , . . .,v M ) (8.4) 

Let us recall that the BA state on the l.h.s. is annihilated by J + for any q (including q 
a root of unity) [3]. Hence eq.(8.4) represents the mixing of two reps with spin J and 
J' = J + 1 into a type I representation. Indeed, applying the mixing rule (8.1) into 
the necessary condition (8.3) for one root at infinity, yields an identity as required: 

1 = g4(J+l) = g 2(J+J'+l) = {gP) 2n = l (g5) 

where n is a suitable (J— dependent) positive integer. 

We can generalize this analysis to any number of vanishing BA z— roots. Let us 
identify the roots going to the origin with z\, Z2, ■ ■ ■ , z r , 1 < r < M. The BAE for 
the remaining non-zero roots z r +\, . . . , zm take the standard form (4.12) valid for a 
BA state formed by M — r B— operators, which has therefore spin J + r. The BAE 
for the vanishing roots take the form 

g 4J+2(r+l) = p TT ^ ~ *k X < J < r (g 6) 

, A I - '-'■■'I 

where the definition of Fj is understood in the limit of vanishing z\, . . . ,z r . We can 
obtain complete agreement with the quantum group mixing rules (8.1) by setting 
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r < p and making the natural choice 

Fj = l, l<j<r (8.7) 

Indeed the presence of r vanishing z— roots imply the mixing of two representations 
with spin J and J' = J + r, so that 

q U+2(r+l) = q 2(J+.r+l) = {q) 2n = 1 (8g) 

with n an integer depending on J and J'. Eq. (8.7) ha a very simple solution for the 
limiting behaviour of the vanishing roots. We find that if 

Zj = uj^ 1 zi , 1 < j < r (8.9) 

with uj r — 1, then identically 

nsin[— 7 + ir(k — ? ) /r 1 
. L r , ^1 = 1, Kj<r 8.10 
( | sm[ 7 + ir(k - j)/r\ > -J- ^ > 

for any value of 7. That eq.(8.10) should hold for generic values of q = e* 7 is necessary 
since we are studying the approach of the roots to the origin when 7/71 tends to a 
rational number. We would like to remark that the limiting behaviour of the ratios of 
the vanishing BA z— roots depend only on their number r and neither on the specific 
rational value of 7 /7T nor on the other non-vanishing roots. 

In the v— plane, the roots v\, i>2, • • • , v r go to infinity as a r— string with spacing 
in/r. For example, when r = 2 we have a pair of complex roots with imaginary 
parts tending to ±27r/4 as their real parts simultaneously diverge. For r = 3 there 
will be a limiting 3-string with a real member and a complex pair at ±«7r/3. For 
these two examples, we have performed an explicit numerical test when 7 — > (ir/3)~ 
and 7 — > (7r/4) ~, respectively, finding perfect agreement with the picture proposed 
here. We shall present more detail on the special cases in the next section. Let 
us anticipate here that the numerical results suggest the following conjecture on the 
behaviour of the singular roots: \zj\ ~ 0(e 1//r ) and Im[t>j — vj(e = 0)] ~ 0(e 1//r ) where 
e = 7 — 7r/ (r + 1). 

When q p = ±1 and the SU(2) q reps divide into type I and type II, it is possible 
to perform the consistent RSOS(p) reduction of the Hilbert space. This consists 
in keeping, among all states annihilated by J + which form the SOS subspace, only 
the type II states. It is known that this corresponds to SOS configurations with 
local heights £ n restricted to the set 1, 2, . . . ,p — 1 [(8.3)6]. In our BA framework, 
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therefore, the RSOS reduction is obtained by retaining only those states with spin 
J < (p — 1)/2, that is with M > (N—p + 1)/2, such that all BA roots z\, Z2, ■ ■ ■ , z r are 
non-zero. This constitues the general and simple prescription to select the RSOSfjo) 
subspace of eigenstates of the evolution operator when q p = ±1. 

We have considered up to here the RSOS reduction at the microscopic level, that 
is for the BA states (4.11) described by the bare BAE (4.12). It should be clear, 
however, that the analysis of the singular BAE solutions when q is a root of unity 
holds equally well for the higher-level BAE (5.18), due to their structural identity with 
the bare BAE. The crucial problem is whether a quantum group reduction carried out 
the higher level would be equivalent to that described above in the "bare" framework. 
This equivalence can actually be established as follows. Suppose that j/tc is irrational, 
but as close as we like to a specific rational value. Then all solutions of the bare BAE 
are regular and can be correctly analyzed, in the limit N — > oo, using the density 
description for the real roots. As described in paper I, this yields the higher-level 
BAE (5.18), those roots are in a precise correspondance with the complex roots of the 
bare BAE. When 7 tends to n/p, then 7 tends to it/ (p — 1) and singular solutions of 
the bare BAE with two or more singular roots are in one-to-one correspondance with 
singular solutions of the higher-level BAE. Thus in this case quantum group reduction 
and renormalization of the BAE commute. The only potentially troublesome cases 
are those of singular solutions with only one real singular root. Indeed the density 
description of the v— roots cannot account, by construction, for real roots at infinity, 
and no sign of the type I nature of the corresponding BA state would show up at 
the higher level. Notice that spin J BA states with a single singular root are mixing 
with spin J + 1 states, so that, by the quantum group mixing rule (8.1), necessarily 
J = p/2 — 1. On the other hand, a direct application of the constraint (8.2) at 
the higher level, that is with the replacement p — > p — p — 1, would overdo the 
job, by incorrectly ruling out all BA states with J = p/2 — 1 = (p — l)/2. To be 
definite, consider N even. Then the type II BA states with J = p/2 — 1, which ar 
e superpositions of SOS configurations with £ n <p— 1 = 2J+1= £n, have an 
effective, higher-level spin J = (i^ — l)/2 = J — 1/2. This higher-level spin does 
satisfy (8.2). The type I states are those in which £ n somewhere exceeds p — 1. In 
particular, the simplest change on a type II configuration, turning it into type I, is to 
flip the oscillating i n close to the right wall (as shown in fig. 6). £n and hence the 
spin J are left unchanged, but clearly £n increases by one, causing J to violate the 
bound (8.2). In other words, RSOS (p)— acceptable BA states with J = p/2 — 1 have 
one kink less than the corresponding SOS states. From the detailed analysis of the 
p = 4 case, to be discussed below, it appears that the removal of one kink corresponds 
to giving infinite rapidity to the hole representing that kink in the higher-level BAE. 

Thank to the possibility of performing the restriction directly at the renormalized 
level, the kink S— matrix of the RSOS(p) model follows by direct restriction on the 
SOS S— matrix given by eq.(7.3). Namely, one must set 7 = n/p, that is 7 = ir/ (p— 1), 
and consider all indices as running from 1 to p — 2. 
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9. The models RSOS(3) and RSOS(4) 



In order to clarify matters about the BA RSOS reduction discussed in the pre- 
vious section, we present here the more details about the two simplest examples, 
when 7 = 7r/3 and 7 = it/4. For these two cases the RSOS reductions correspond, 
respectively, to a trivial one-state model and to the Ising model (at zero external 
field and non-critical temperature). Let us recall that the light-cone approach yields 
in the continuum limit massive field theories with the same internal symmetry of the 
corresponding critical regimes. Therefore the RSOS (4) reduction of the light-cone 
6V model coincides with the preserving perturbation of the c = 1/2 minimal 

model (which is obtained upon quantum group restriction of the critical 6V model). 

9.1. The case p = 3 

From the RSOS viewpoint, the case 7 = tt/3 is particulary simple. Since p = 3 the 
local hight variables £ n can assume only the values 1 and 2. Then the SOS adjacency 
rule \£ — £'\ = 1 implies that the global configuration is completely determined once 
the height of any given site is chosen. In our light-cone formulation the first height 
on the left is frozen to the value 1 (see eq.(4.18) and fig. 2), so that the restricted 
Hilbert space contains only one state: an SU (2) q singlet when N is even and the spin 
up component of a J = 1/2 doublet when N is odd. 

For 7 = 7r/3 the BAE still have many solutions which reproduce the full SOS 
Hilbert space. Indeed, from the SOS point of view, nothing particular happens 
when 7 — > 7r/3. However, at this precise value of the anisptropy, only one BAE 
solution corresponds to the unique type II representation: for N even (odd) it con- 
tains N/2 (N/2 — 1/2) non-zero roots in the z— plane. All other BAE solutions with 
M = [N/2] contain at least one vanishing z— root and correspond to type I SOS 
states. Moreover, there unique RSOS state is the ground state of the SOS model and 
is formed by real positive roots labelled by consecutive quantum integers (no holes). 
We reach therefore the following rather non-trivial conclusion: the complicated sys- 
tem of algebraic equations (4.12) admit, for g 3 = ±1 and w real, one and only one 
solution with M = [N/2] non-zero roots within the unit circle. In addition, these 
roots are real and positive. For few specific choices of N we also verified this picture 
numerically. 

9.2. The case p = 4 

The local height variables can assume now the three values £ = 1,2, 3. However, 
each configuration can be decomposed into two sub-configurations laying on the two 
sublattices formed by even and odd faces, respectively. On one of the two, the SOS 
adjacency rule \£ — £'\ — 1 freezes the local heights to take the constant value £ = 2. 
Then on the other sub-lattice we are left with two possibilities, £ = 1 and 3. Moreover, 
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the interaction round-a-site of the original model reduces in this way to a nearest- 
neighbor interaction in the vertical and horizontal directions. The framework is that 
of the Ising model. 

To obtain the standard Ising formulation, we can set 

a n (t)=£ 2n (t)-2 (9.1) 

where the numbering of the lattice faces can be read from fig. 2. The fixed b.c. on 
the £ n now correspond to 

ao(t) = -l, a R (t) = (-l) J - 1 (N = 2R) 
<7 (*) = -l, a R (t) = ±l (N = 2R+1) 

where J = 0, 1 for iV even, due to the bound (8.2). For iV odd we must consider only 
the possibility J = 1/2, which implies tjq — 2, since the line of half-plaquettes on 
the extreme right belong to the frozen sublattice. Then £n-i is left f ree to fluctuate 
between 1 and 3, leading to free b.c. on a. The matrix elements of the unit time 
evolution operator U can be calculated from the explicit form of the SOS weights 
(4.27)-(4.28). They can be written in the "lagrangian" form 

m (a' , a[, . . . a' R \ U \a , a u . . . a R \ = (9.3) 

where 

R-l R-l 

L{t) =(3 V J2 Vn(t)a n (t + l)-f3 h J2 Vn(t)a n+1 (t) + const (9.4) 

n=l n=0 

and 

V = J(tt + 2ilntanh26) , (3 h = arctantanh29 (9.5) 

Alternatively, standard simple manipulations allow to rewrite U explicitly in terms 
of Pauli matrices 

jj =e -i/3hH2 e -iPhH! 

^ ?4 (9.6) 

n—1 n—Q 

In either cases, one sees that the complex Boltzmann weights formally belong to the 
critical line 

sin 2f3 v sin 2(3 h = 1 (9.7) 

Of course, this follows from the original definition of the light-cone 6V model in terms 
of complex trigonometric Boltzmann weights, which, under the replacement O — > iO, 



43 



would correspond to the critical standard 6V model. Nevertheless, just as for the 
vertex model, also in this "light-cone" Ising model, a massive field theory can be 
constructed in the Re 6 — > +00 limit. 

In the BA diagonalization of U, we must restrict ourselves to the BAE solutions 
with M = R or M = R - 1 for even JV and M = R for odd JV. For sufficently 
large even JV, the ground state has R real roots with consecutive quantum integers. 
Actually, as already stated above, this is a general fact valid for all RSOS(p) models. 
Namely, the infinite volume ground state of th SG model, of the SOS model and of all 
its restrictions RSOS(p) is the same f.b.c. BA state. It is the unique SU(2) q singlet 
with all real positive roots and no holes. It is described in the thermodynamic limit 
(JV — * 00, a fixed) by the density of roots given in eq.(5.7) of paper I. 

Excited states with J = have an even number of holes and a certain number of 
complex roots such that no z— root lays at the origin. For instance a two-particle state 
contains two holes (holes are naturally identified with the particles) and a two-string 
with imaginary parts ± [71-/8+ corrections exponentially small in JV] laying between 
the two holes. This state (that is this precise choice of quantum integers) is just the 
two-particle state of the SG or SOS model, with 7 fixed to the precise value 7r/4. 
We have explicitly checked, by numerically solving the BAE for various values of JV, 
that indeed all z— roots stay away from the origin as 7 crosses 7r/4 while the quantum 
integers are kept fixed to the two-hole configuration. Now consider a state with four 
holes. Apart from the multiplicity of the rapidity phase space, there are two distinct 
type of such states: in v— space one contains two two-strings, while the other contains 
one sigle wide pair, that is a complex pair with imaginary part larger than 7. This 
situation holds for generic values of ^/ir and simply reflects the fact that the holes are 
SU(2) q doublets. The crucial point is that, as 7 reaches it/ 4 (from below), the wide 
pair moves towards infinity, while its real part gets closer and closer to the (diverging) 
value of the largest real root and the imaginary part approaches it/3. This picture is 
confirmed by a careful numerical study of the BAE and is in perfect agreement with 
the general picture presented in sec. 5. In addition, there exists numerical evidence 
that the real parts diverge like log(7r/4 — 7)~ 1//6 while the imaginary part of the wide 
pair goes to 7r/3 like (tt/4 — 7) 1 / 3 . When 7 > 7r/4, the v — root corresponding to the 
largest quantum integer has the largest real part, but is no longer real, having an 
imaginary part equal to it/ 2. As 7 — > (tt/4) + , this real part again diverges together 
with real part of the wide pair. Hence the four-hole state with a wide pair is type I 
when 7 = 7r/4. On the other hand, the four-hole state with two two-string is type 
II, since all its v— roots stay finite. 

From the study of the two- and four-hole states, we are led to the following 
general conjecture: in the BA framework, the RSOS(4) J = states are all and 
only the states with v — 2k holes and k two-strings. Then, in the higher-level BAE 
(5.18), we recognize this state as that corresponding to the unique solution with k real 
X~ roots. In other words, this BA state is completely determined once the location 
of the holes (that is the rapidities of the physical particles) is given. 
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Consider now the states with J = 1. In this sector, the lowest energy type II state 
contains exactly one hole. This is an unusual situation for BA systems on lattices 
with even N, where holes are always treated in pairs. As long as 7 7^ it/ 4, the same 
is true in our f.b.c. BA: the lowest energy state with J = 1 contains two holes, in 
agreement with the interpretation of holes as SG solitons with quantum spin 1/2. For 
7 < 7r/4 all roots are real. For 7 > 7r/4 the root vjv/2-1 corresponding to the largest 
quantum integer In/2-1 = N/2 + 1 aquires an imaginary part n/2 (see the appendix 
for details). But when 7 = 7r/4 then Revjv/2-1 — +00, and the J = 1 two-hole state 
is mixed with some J = 2 state into a type I representation. It does not belong to 
the RSOS(4) Hilbert space. This picture can be easily verified numerically. 

To prevent the largest root v N/2-1 from diverging, it is sufncent to consider a 
J = 1 state with only one hole and ijv/2-i = N/2. From the point of view of the SG 
or SOS models this one-hole state has a cutoff-dependent energy which diverges as 
a -1 in the continuum limit. If x\ is the position of the single hole, the energy relative 
to the ground state reads 

E - E = a _1 e(xi) + na^ 1 (9.8) 

where the renormalized energy function e(x) is given in eq.(6.1). Notice that this 
result coincides with the limit xi — > 00 of a two-hole state (eft. eq.(6.1)). Strictly 
speaking, however, the density method leading to (6.1) has no justification when "one 
hole is at infinity". We trust eq.(9.8) nonetheless because it passes all our numerical 
checks. In the continuum limit a — * 0, — » 00, 4a _1 e~ 7r<9//7 = m (fixed), this one- 
hole state is removed from the physical spectrum, as required from the SG point of 
view: holes are spin 1/2 solitons while here J = 1. 

The preceding discussion easily extends to generic states in the J = 1 sector 
containing an odd number of holes and a given set of complex pairs. The infinite- 
volume energy of these multiparticle states is given by 

V 

E = EQ + a~ l ^e(x h ) + na" 1 (9.9) 
h=l 

with the same divergent constant 7ra _1 appearing irrespective of the physical content 
of the state. Our conjecture for the J = states naturally extends to these J = 1 
states: if there are v = 2k + 1 holes, then the v— roots are all finite, implying that 
the state is type II, provided there are also k two-strings. Then the higher-level BAE 
imply that these states are completely defined by the hole rapidities. To retain these 
type II J = 1 BA states in the RSOS(4) model, an extra J— dependent subtraction 
is necessary to get rid of the divergent constant. Namely, for even and J = or 1, 
we set 

Hrsos(4) = Hsosil = tt/4) - tto" 1 J (9.10) 
where H SO S is given by eqs. (4.10), (4.14), (4.22)-(4.28). 
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Altogether, we see that the BA picture for the excitations of the RSOS(4) model 
is fully consistent with the kink interpretation for the holes. Indeed J = corresponds 
to even Dirichlet b.c. for the Ising field, while J = 1 corresponds to odd Dirichlet 
b.c. (eft. eqs. (9.2)). In the Ising model the kink description is valid below the 
critical temperature, with the disorder field as natural interpolator for the kinks. In 
this case the S— matrix must be —1 (it is +1 when the asymptotic particles are the 
free massive Majorana fermions) and this is exactly what follows from eq.(7.3) upon 
setting 7 = ir/3 and restricting all indices to run from 1 to p — 2 = 2. 

The analisys of the BA spectrum for odd iV does not contain real new features. 
Now J is fixed to 1/2 and the b.c. on the Ising field are of mixed fixed-free type, 
as shown in the second of eqs. (9.2). The lowest energy state corresponds to the 
BAE solution formed by real roots and one single hole as close as possible to the 
the v— origin. In other words, the quantum integers are given by Ij — j + 1, for 
j — 1, 2, . . . , (N — 1)/2. By letting this hole to move away along the positive real axis, 
we reconstruct the energy spectrum of a one-particle state. Of course, to compare 
this state to the global ground state, which contains no holes, a judicious choice of 
the odd value of N is required. If in the ground state N = 2R, then the new state is 
indeed an excited state if we choose N = 2R — 1, rather than N = 2R+ 1, because 
of the antiferromagnetic nature of the interaction. 
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In this appendix, we present for completeness a (rather non-standard) treatment 
of the BAE (4.12). As explained in paper I, it is convenient to first rewrite them in 
the p. b.c. form 



sinh(Aj - + ry/2) sinh(Aj + 6 + ry/2) 



i N 



sinh(Aj - - ry/2) sinh(A J - + 6 - ry/2) 

sinh(Aj — Afc + ij) 



sinh(2Aj + ry) 
sinh(2Aj — ry) 



n 



k=-M+l 



sinh(Aj — Afc — ry) 



(A.I) 



where the numbers A_m+i, ■ ■ ■ , are related to the v— roots by 



A i = v i = 



-A 



j = l,2,...,M 



(A.2) 



and Revj > 0, thanks to the symmetries of the BAE (4.12). Let us concentrate our 
attention on those BAE solutions which are mostly real, that is those which contain 
an arbitrary but fixed number of complex pairs interspersed in a sea of order N of real 
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roots. The reason for this restriction will become clear in the sequel. The counting 
function associated to a given solution v i, . . . , vm is defined to be 

Z N (x;v\,..., v M ) =0 7 /2(^ + @) + 7 /2(^ - B) + N~ 1 (p 1 (2x) 



j=-M+l 



where 



, , . .. sinhfia + x) ,. A . 

Mx) = % log — - y. '- A.4 

smh(m — x) 

The logarithmic branch in eq.(A.4) is chosen such that <f> a (x), and as direct conse- 
quence Zn(x ; vi, . . . , vm), are odd. 

The BAE (A.l) can now be written in compact form 

Z N (X j ;v 1 ,...,v M ) = 2ttN- 1 I j , j — —M + 1, —M — 2,...,M (A.5) 

where the quantum integers Ij entirely fix the specific BAE solution, and therefore 
the BA eigenstate, and by construction satisfy = —Ij for j — 1,2, ... , M. Notice 
also that by definition Zj^(0 ;v\, . . . , vm)— 0, but A = is not a root, due to (A. 2). 
Rather, A = is always a hole, from the p.b.c. point of view. 

To begin, consider the case when iV is even, M = N/2 (i.e. J = 0) and all roots 
are real. This state (the groud state for even N) is unambiguously identified by the 
quantum integers 

Ij=3 j = l,2,...,N/2 (A.6) 

and the corresponding counting function is indeed monotonically increasing on the 
real axis, justifying its name. To our knowledge, the existence itself of this BAE 
solution has not been proven in a rigorous analytical way. But it is very easy to 
obtain it numerically for values of iV in the thousands and precisions of order 1CP 15 
on any common workstation. 

Next consider removing J roots from the ground state. J is the quantum spin of 
the corresponding new BA state. For 7 suflicently small, one now finds that 

N/2 + J < Z N (+oo;v h ...,v M ) < N/2 + J + 1 (A.7) 

so that, together with the actual roots satisfying eq.(A.5), there must exist positive 
numbers x\, X2, • • • , x v , with v > 2J , satisfying 

Z N (x h ;v 1 ,...,v M ) = '2nN~ 1 I h , h = l,2,...,v (A.8) 

where the 1^ are positive integers. If Zj^ix) is monotonically increasing, then nec- 
essarily v = 2 J and the integers . . . , ij\r/2-j> h, ■ ■ ■ , hj} are all distinct. The 



47 



numbers x\, X2, ■ ■ ■ , x v are naturally called holes. For J held fixed as N becomes 
larger and larger, it is natural to expect that the counting function is indeed mono- 
tonically increasing, and numerical calculations confirm this expectation. For larger 
values of 7 the situation becomes more involved. Numerical studies show that, first 
of all, Z n (x) develops a local maximum beyond the largest root, while still satisfying 
the bounds (A. 7), as 7 exceeds a certain (J— dependent) value. For even larger values 
of 7, the asymptotic value of Z^(x) becomes smaller than Z* = 2n(l/2 + J/N), but 
its maximum stays larger than Z*, provided there is indeed a root corresponding to 
N/2 + J (i.e. Im = N/2 + J). Up to now, vm is obviously located where Z^(x) 
reaches N/2 + J from below, and one could say that there exists an extra hole X2J+1 
further beyond, where Zn(x) reaches Z* from above. When 7 reaches a certain crit- 
ical value, the local maximum lowers till Z*. At this point the root and the extra 
hole exchange their places, and for sligtly larger values of 7 the hole is located where 
Zn(x) reaches Z* from below, while the root lays further away, where Zjy(x) reaches 
Z* from above. As the numerical calculations show, however, this extra hole with 
the same quantum integer Im = N/2 + J of the largest root hole is spurious, since no 
energy increase is really associated to its presence. When none of N/2 — J root has 
N/2 + J as quantum integer, then Zn(x) does never reach Z* for sufficently large 7, 
and, strictly speaking there are only 2 J — 1 holes. This time, however, we find that 
the energy increases with respect to the ground state in the same way as if there was 
"a hole at infinity", that is a hole beyond the largest root. Thus v can always be 
regarded to be even, when N is even. 

As an important example consider the definite choice J = 1. Then for 7 < 7r/6 
we find Z^(+oo) > Z*, and there are two holes, with 1 < I\ < I2 < N/2 + 1. Assume 
I2 < N/2 + 1 and consider the interval 7r/6<7<7r/4. The counting function has a 
maximum Z max (larger than Z*) situated to the right of the largest root f/y/1-1, as 
long as 7 is smaller than the critical value 7* at which the maximum lowers to Z*. 
For 7 > 7*, the maximum is still larger than Z* but is located to the left of ujv/i-i- 
In any case, for n/6 < 7 < 7r/4 the asymptotic value Z/y(+oo) is smaller than Z*. 
When 7 — > (7r/4) ~, then the largest root as well as the maximum Z max are pushed to 
infinity, and Zjq(x) is once again monotonic with Zjv(+oo) — Z* . As 7 exceeds 7r/4, 
the last root v jv/1-1 passes, through the point at infinity, from the real line to the line 
with Imu = 7r/2. This pictures generalizes to arbitrary J with the two special values 
7 = 7r/6 and 7 = n/4 replaced, respectively, by 7 = re/ (4J + 2) and 7 = ir/ (2 J + 2). 
In fig. 7 the salient portion of the numerically calculated counting function is plotted 
for J = 1, iV = 64 and a specific choice of Ii, J2. In this case we approximatively 
find 7* ~ 0.2br. 

Finally, consider a BAE solution containing, in addition to a number of order N 
of real roots, also a certain configuration of complex roots. In the v— space, these 
complex roots appear either in complex conjugate pairs or with fixed imaginary part 
equal to in/2, so that the counting function is real analytic: Z^(x) = Z^(x). More- 
over, it is fairly easy to show, by looking at the value of the counting function a real 
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infinity, that the presence of complex roots implies the existence of holes in the sea 
of real roots. For our next porposes, we shall now consider 7/7? irrational, so that 
all v— roots are finite. Denoting with u q , q = 1,2, ... , M c the values of the complex 
roots and with M r = M — M c the number of real roots, we now write the derivative 
of the counting function as 

2np N (x) = Z' N (x) = F^(x) + N- 1 [Fftx) + F' c {x) - F' h {x)] + (K * p 5 )(x) (A.9) 

where 

Fq(x) =0 7/2 (x + 0) + 7/2 (a; - 0) 
F (x) =0 7 (2x) 

F c( x ) = ^2 [<i>i( x ~ u q) + <l>i( x + u q)] 

q=1 v (A.10) 

F h {x) =(j)j(x) + ^ [^>l{ x ~ X n) + 07^ + x n)\ 
n=l 

M r V 

Nps(x) = y^ j 6(x- Xj) + 7 (x) + ^2 [H x - x h) +d( x + x h)\ 

j=l n=l 

and K* is the convolution defined by 

+00 

(K*f)(x)= J dy^(x-y)f(y) (A.11) 

— 00 

The so-called "density approach" consists in replacing, as A" — > 00, both p^ and pg 
with the same smooth function p, representing the density of roots and holes on the 
real line. This function is therefore the unique solution of the linear integral equation 
(eft. eq.(A.9)) 

2np = Fo + N~\F + F c -F h ) +K*p (A.12) 

which can be easily solved by Fourier transformation. Combining this equation with 
eq.(A.9), we now obtain, after some simple manipulations 

PN = P + G * (pn - Ps) (A.13) 

where G* = (2% + A') -1 * K* stands for the convolution with kernel 

G ( x ) = f dk c *kx sinh(7r/2- 7 )A: 
v ; J 2tt sinh(7r - 7 )A;/2 cosh 7 A;/2 v ' ; 

Finally, a simple application of the residue theorem to the analytic function pjv(l — 
e -iNZ N ^-i pj ug an i n t e g ra tion by parts lead to the following formal nonlinear integral 
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equation for the counting function 



Z}\[ — Z + G * Ljv 



(A.15) 



where 



L N (x) = -iN' 1 log 



]_ _ e iiVZiv(a;+iO) 



(A.16) 



X _ e -iNZ N (x-iO) 



and Z is the odd primitive of 27rp, namely 



^ + 7V- 1 (2tt + X) * (F + F c - F h ) 



(A.17) 



with 



^(x) = ((2?r + X) * Fq) (x) = 2 arctan 



cosh7r6/7 



sinh7rx/7 



) 



(A.18) 



Eq.(A.15) is a formal integral equation since the knowledge of the exact position of 
holes and complex roots is required in Z. It becomes a true integral equation for the 
ground state counting function. In any case it is an exact expression satisfied by Zn 
where the number of site N enters only in an explicit, parametric way, except for the 
positions of the holes and of the complex roots. After having fixed the corresponding 
quantum integers, these parameters retain an implicit, mild dependence on N, for 
large N. 

We shall now show that eq.(A.15) is very effective for establishing the result (6.6), 
which is of crucial importance for the calculation of the S— matrix. It is sufncent to 
check that the second nonlinear term in the r.h.s. of eq.(A.15) is indeed of higher 
order in iV _1 relative to the first. To this purpose observe that the integration contour 
of the convolution in eq.(A.15)can be deformed away from the upper and lower edges 
of the real axis, since for sufficently large N no complex roots can appear in the whole 
strip |Imx| < 7/2 and G(x) is analytic there. Indeed Z^ tends to Z^ as N — > 00, 
and one can eplicitly check that the imaginary part of Z^ is positive definite for 
< Imx < 7/2 and negative definite for > Imrr > —7/2. This also implies that 
the contribution to the convolution integral is exponentially small in N for all values 
of the integration variable (let's call it y) where ImZoo of order 1. For |Rey| of 
order log A, we find Ym. Z^y) of order A -1 , so that the nonlinearity L^v, rather than 
exponentially small, is also of order A -1 . But now the exponential damping in y of 
the kernel G(x — y) guarantees that the convolution integral is globally of order iV -2 
or smaller, provided \x\ is kept smaller than (7/77) logiV. Hence we can write 



Finally, the coefncent of the iV 1 term of Z, in eq.(A.17), can be calculated in the 



Z N = Z + 0(N- 2 ) 



(A.19) 
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N — > oo limit, with the techniques described at length in paper I. After some straight- 
forward albeit cumbersome algebra, this yields 

Jim ((2tt + K) * (F + F c - F h )) (x) = F(x ; x u . . . , x M - Xl , . . . , qx M ) (A.20) 

N^oo 

where the higher-level quantity F is defined in eq.(6.8). Together with eqs. (A. 19) 
and (A. 17), this proves eq.(6.6) of section 3, as claimed. 

Let us now consider the problem of calculating the energy of a given BA state, 
through eq.(4.15). We rewrite first the "bare" energy function Qq{x) as 

eo(a;) = -2tt + <f)^ 2 (x + 0) - 7/2 (:r - 0) (A.21) 

Then we calculate 

3=1 -oo 1 =1 ^=1 

+?° M c 
= \N J p(x)0 7/2 (x + 0) + ^ 7/2 (« g ) -J2^/2(xh) ~ 50 7 /2(O) 



-oo 



+00 



q=l h=l 



+±iV J \ps(x) - pn(x)] ((1 -G)* 7 / 2 ) (s) 



(A.22) 

Through the residue theorem and an integration by parts, the last term can be trans- 
formed, as done before for the counting function, into an integral of the nonlinear 
term L/y, namely the integral 



+oo 



2 iv 



L N (x) 



7 cosh nx/'y 



—oo 



By the same argument used above, this last expression is globally of order AT -1 . 
Finally, inserting the explicit form of the continuum density p(x) into eq.(A.22) and 
recalling eqs. (4.15) and (A.21), for the energy we obtain eq.(6.1) of the main text 



E = Eo + a^J^^^+Oia^N- 1 ) (A.23) 
h=l 
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where 



E = - na 1 N + 



oc 




dk sinh(7r — •y)k/2 sin kO 



k sinh7r/c/2 cosh -fk/2 



2Ncosk& + 1 + 



sinh(7T - 3 7 )A;/41 (A.24) 
sinh(7r — 7)fc/4 



o 



is the energy of the ground state. 

We would like to close this appendix with a comment on the limitations of the 
density approach, where one deals only with the solution p of the linear equation 
(A. 12). Regarding p(x) as the actual density of real roots and holes in the N — > oo 
limit, it is natural to use it to replace summations with integrals. What one learns 
from the exact treatment presented above as well as from computer calculations, is 
that the error made in such a replacement depends crucially on the large x behaviour 
of the quantity which is to be summed. This error is down by N^ 1 only when there 
is exponential damping in x. This means, for instance, that the integral of p(x) does 
not reproduce in general the exact number of real roots and holes, but rather some 
7— dependent quantity close to it. Misundertanding this for the actual number of real 
roots and holes would lead to the absurd result that the holes have a 7— dependent 
value of the SU(2) q spin, which is instead necessarily integer or half-integer and, in 
the particular case of the holes, just 1/2 for any value of 7. 
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11. Figure Captions 
Fig.l. Graphical representation of the i?-matrix. 

Fig. 2. The standard row-to-row monodromy matrix. The numbers from 1 to N label 
the vertical spaces and each vertical line represents a couple of free indices in 
the corresponding vertical space. Indices on the internal horizontal lines are 
summed over. 

Fig.?>. The doubled monodromy matrix used to describe systems with fixed boundary 
conditions. The correct contraction over the internal indices is dictated by the 
position of the arcs. 

Fig A. Graphical representation of the compatibility relations (2.12) between the two- 
body scattering and reflections on the left wall. 

Fig. 5. Graphical representation of the compatibility relations (2.15) between the two- 
body scattering and reflections on the right wall. 

Fig. 6. The unit time evolution operator U(Q) for even N (top) and odd N (bottom). 
As usual, indices over internal lines are summed over. 

Fig. 7. Reconstruction of the diagonal light-cone lattice through powers of the unit 
time evolution operator (case N even). 

Fig. 8. Assignement of the local height variables on the plaquettes cut by the time- zero 
and time-one lines. As is evident from the figure, these two lines sandwich the 
unit time evolution operator [/(©). 

Fig.9. The Bratteli diagram corresponding to the case N=6. Right-moving paths on 
this diagram arriving at height j define the basis in the space of irreps of weight 
3- 

Fig.10. The "ground state dominating" configuration of local height variables. 

Fig. 11. Graphical representation of the vertex-face correspondence. The strip sand- 
wiched between the two time lines represents the unit time evolution. 
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Fig. 12. Graphical representation of the one-soliton operator corresponding to the par- 
ticle with rapidity 62 in a system with v — 4 particles. At each intersection, 
the appropriate two-body S— matrix acts. In the collisions against the wall, 
62 is flipped and the two-by-two matrix (72 (±02) ac ts on the internal states of 
the soliton. The choice of 62 as largest rapidity is done purely for graphical 
convenience. 

Fig. 13. One of the local height configurations that dominate the one-hole BA state, of 
spin J = 1. 

Fig. 14. The "renormalized" version of Fig. 13. 

Fig. 15. Flipping the last portion of the configuration from the solid to the dotted line 
does not change the quantum spin J = p/2 — 1, but tranform the type II state 
into type I. 

Fig. 16. Plot of the N/2n times the counting function Zjq{x) versus tanhx for N = 64, 
J = 1, = .15 and various values of 7 from n/Q to 0.9997r/4. The quantum 
integers associates to the two holes are I\ — 8 and I\ = 21. The critical value 
of (N/2ir)Z N is N/2 + J = 33, while that of 7 is, roughly, 7* = 0.21vr. 
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